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The  objective  of  this  research  program  has  been  to  investigate  the  potential  of 
new  fiber  optic  configurations  for  inertial  rotation  sensing  so  as  to  extend  the 
sensitivity  and  stability  of  such  systems.  In  particular,  the  research  has  been 
aimed  at  the  design  and  construction  of  reentrant  Sagnac  systems  permitting  re¬ 
circulation  of  the  optical  signals. 


There  have  been  two  main  activities  in  the  program,  (1)  the  design  and  partia 
construction  of  an  all-fiber  Sagnac  rotation  sensor  eliminating  all  previously — [ 
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-used  bulk  optical  components,  lenses,  mirrors,  prisms,  etc.  Bulk  components  do 
not  have  sufficient  precision  and  impose  alignment  problems  so  as  to  limit  the 
sensitivity  and  recirculation  capability  of  Sagnac  sensors.  The  all  fiber  de¬ 
sign  was  made  possible  by  the  development  in  this  laboratory  partially  under 
other  contract  support  of  a  set  of  all  fiber  components,  directional  couplers, 
polarizers,  and  polarization  controllers,  which  have  been  incorporated  in  the 
all  fiber  design.  This  configuration  should  provide  greatly  increased  sensiti¬ 
vity  and  device  compactness  and  ruggedness.  The  second  set  of  activities  in  the 
program  has  been  to  investigate  various  limitations  on  the  sensitivity  (minimum 
detectable  rotationrate)  of  Sagnac  systems.  One  such  limitation  is  the  extrane¬ 
ous  noise  known  to  exist  in  all  Sagnac  systems.  We/have  found  particularly, 
that  the  effects  of  Rayleigh  backscattering  from  fiber  inhomogeneities  of  phase 
detection  in  Sagnac  systems,  as  one  of  the  principle  limitations  on  sensitivity. 
Investigation  of  various  possible  means  of  circumventing  this  very  fundamental 
noise  source  are  underway.  The  second  major  effect  on  Sagnac  sensitivity  comes 
from  signal  polarization  and  the  effects  of  fiber  properties  on  polarization 
conversion  along  the  fiber  path.  A  detailed  mathematical  formalism  has  been  de¬ 
veloped  to  treat  the  effects  of  polarization  conversion  on  ultimate  Sagnac  sen¬ 
sitivity.  As  a  result  of  this  mathematical  treatment,  a  particular  form  of 
source  polarization  and  a  specific  set  of  operating  conditions  for  a  Sagnac  sen¬ 
sor  has  been  suggested  which  should  minimize  the  deleterious  effects  of  polari¬ 
zation  conversion.  This  configuration  has  been  tested  and  has  worked  quite 
successful ly. 
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I.  INTRODUCTION 


A.  OBJECTIVE 

The  objective  of  this  research  program  has  been  to  investigate  the  po¬ 
tential  of  new  fiber  optic  configurations  for  inertial  rotation  sensing, 
so  as  to  extend  the  sensitivity  and  basic  stability  of  such  systems.  The 
approach  has  involved  basic  device  research,  including  both  theoretical 
analysis  and  experimental  demonstrations  for  the  purpose  of  improving 
future  generations  of  operational  rotation  sensing  systems.  We  have  intro¬ 
duced  new  concepts  in  the  use  of  optical  fibers  in  rotation  sensing  loops, 
specifically  the  recirculation  of  optical  pulses  in  these  loops  (reentrant 
Sagnac  system),  and  eventually  active  operation  of  these  loops  (i.e.,  the 
incorporation  of  an  amplifier  in  the  loop),  aimed  at  increasing  the  sensi¬ 
tivity,  stability,  and  optical  integration  times  of  these  devices  by  very 
large  factors. 

B.  PROGRAM  SUMMARY 

This  program  has  had  two  main  sets  of  activities  over  the  past  year. 

The  first  was  concerned  with  continuation  of  the  previous  work  under  the 
AFOSR,  aimed  at  the  design  and  construction  of  an  active  reentrant  Sagnac 
(ARS)  pulsed  all-fiber  optic  system  for  inertial  rotation  sensing.  This 
system  is  described  in  Section  II.  The  current  phase  of  the  work  is  lead¬ 
ing  to  a  passive  all-fiber  reentrant  system  (PRS).  There  are  no  all-fiber 
reentrant  gyros  of  any  kind  in  existence  at  this  time,  to  our  knowledge. 

The  second  activity  has  been  concerned  with  examining  various  problems 
related  to  the  sensitivity  (minimum 
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systems.  R  &  D  on  devices  in  which  CW  light  waves  make  a  single  pass 
round  a  multi -turn  fiber  optic  sensing  loop  is  being  carried  on  in  a 
number  of  laboratories.  At  the  start  of  the  present  reporting  period, 
all  such  devices  were  limited  to  sensitivities  far  below  theoretical 
predictions  and  there  v/as  an  urgent  need  to  overcome  these  noise  bar¬ 
riers.  This  problem  became  an  important  part  of  our  activity.  In  this, 
we  are  largely  concerned  with  studying  the  effects  of  various  sources  of 
noise  in  the  fiber  and  the  effects  of  birefringence  and  nonreciprocity 
in  the  fiber. 

Our  activity  on  these  two  sets  of  problems  is  described  briefly  be¬ 
low  and  in  greater  detail  in  the  body  of  this  report  (Section  II  on  Re¬ 
circulating  Sagnac  System  and  Section  III  on  New  Techniques  for  Fiber 
Gyroscopes) . 

Obviously  methods  for  increasing  the  sensitivity  and  accuracy  in  the 
measurement  of  extremely  small  rotation  rates- are  of  great  importance  to 
Air  Force  objectives  in  the  fields  of  inertial  navigation  and  geodesy. 

It  is  also  relevant  to  point  out  that  fiber  sensors  are  growing  increas¬ 
ingly  important  for  measuring  a  variety  of  environmental  effects  including 
strain,  acoustic  waves,  current,  magnetic  field,  temperature  and  so  forth 
as  well  as  inertial  rotation.  Typically,  such  devices  require  fiber  inter¬ 
ferometers  for  measuring  small  optical  phase  shifts,  and  these  interfero¬ 
meters  all  require  single  mode  fibers.  Consequently,  any  research  on  fi¬ 
ber  rotation  sensors  in  our  present  program  and  the  developing  of  components 
will  lead  to  new  techniques  for  single  mode  fiber  systems  which  are  important 
for  these  other  sensors.  The  results  also  may  be  of  importance  for  signal 
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C.  REVIEW  OF  CURRENT  ACTIVITIES 


The  ultimate  objective  of  the  program  is  an  active  system  incorpora¬ 
ting  an  amplifier  in  the  rotation  sensing  loop.  Our  activity  up  to  now, 
however,  has  involved  rotation  sensing  by  passive  operation  of  pulsed  re¬ 
circulating  fiber  optic  loops,  or  passive  reentrant  Sagnac  systems  (PRS), 
that  is,  with  no  amplifier.  Using  standard  optical  bench  components  it 
has  been  possible,  prior  to  this  reporting  period,  to  observe  five  recir¬ 
culation  pulses  in  the  fiber  loop.  Such  operation  multiplies  the  effective 
length  of  the  fiber  by  the  number  of  recirculations  and  increased  the  rota¬ 
tion  sensitivity  of  the  system  in  accord  with  the  theory.  All  this  refers 
to  PRS  operation  as  described  above  where  the  system  is  operated  in  a  re¬ 
circulating  mode  but  passively.  These  results  showed  that  standard  optical 
bench  components  do  not  have  sufficient  precision  to  align  the  system  for 
simultaneous  optimum  pulse  recirculation  and  rotation  sensitivity.  It  is 
possible  that  more  precise  optical  components  could  overcome  these  problems 
and  improve  performance.  However,  our  experience  indicates  the  required 
precision  and  quality  would  be  difficult  to  achieve.  More  important,  how¬ 
ever,  we  have  developed  a  much  superior  and  simpler  approach  which  will  be 
described  below. 

In  principle  the  mode  and  alignment  problems  with  bulk  recirculating 
components  could  be  completely  circumvented  by  the  use  of  integrated  {single 
mode)  fiber  components.  No  such  components  with  sufficiently  high  perfor¬ 
mance  and  low  enough  loss  were  available  when  the  program  was  started.  Sub¬ 
sequently,  however,  under  a  companion  program  (JSEP)  in  the  laboratory  con¬ 
cerned  with  optical  signal  processing,  we  developed  single  mode  fiber  optic 
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directional  couplers  which  appear  to  satisfy  all  the  requirements  for  a 
PRS.  In  addition,  other  fiber  components,  polarizers,  polarizer  control¬ 
lers  and  so  on  have  been  also  developed  under  another  contract  and  these 
eliminate  all  the  mode  and  alignment  problems  encountered  with  hybrid 
systems.  They  provide  miniaturization,  ruggedness  and  decreased  system 
expense.  We  have  designed  and  tested  various  components  and  sub-systems 
to  be  incorporated  into  an  all -fiber  recirculating  system.  Because  of  the 
unavailability  of  some  particular  optimum  fibers  we  were  not  able  to  com¬ 
plete  the  system,  but  the  design  is  completed  and  the  fabrication  v/as  well 
along  by  the  end  of  the  reporting  period.  These  various  components,  the 
design  of  the  system  and  expected  performance  are  described  in  Section  II 
below. 

It  is  perhaps  useful  to  mention  here  that  such  components  have  been 
incorporated  in  a  CW  nonrecirculating  Sagnac  system.  The  system  and  the 
components  have  operated  very  satisfactorily  and  have  given  higher  sensiti¬ 
vity  in  rotation  sensing  than  anything  else  that  has  been  achieved  in  a 
fiber  system.  It  must  be  emphasized,  again  however,  that  this  is  not  a 
recirculating  system  which  limits  its  ultimate  sensitivity.  This  was  done 
under  the  auspices  of  another  contract,  but  has  provided  very  useful  design 
information  and  operating  experience  which  is  very  relevant  to  the  pulsed 
recirculating  system  which  is  now  under  construction  under  this  program. 
Conversely  many  of  the  design  features  and  operating  characteristics  of 
this  CW  Sagnac  system  are  based  on  data  and  insights  gained  from  other 
AFOSR  work  on  recirculating  systems.  The  existence  of  both  programs  in 
the  laboratory  with  different  objectives  has  been  of  great  mutual  benefit. 
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As  to  the  second  set  of  activities  referred  to  above,  which  are  con¬ 
cerned  with  sensitivity,  basically  we  have  worked  in  two  main  areas.  One 
dealt  with  a  previously  unsuspected  major  source  of  noise  in  the  systems. 
This  is  Rayleigh  scattering  from  inhomogeneities  in  the  fiber  which  produces 
a  noise  level  much  higher  than  the  previously  estimated  theoretical  lower 
limits.  This  work  is  discussed  in  further  detail  in  Section  III,  including 
various  approaches  we  have  developed  to  reduce  this  noise  which  we  are  now 
investigating  under  the  continuation  of  the  contract.  The  proposed  experi¬ 
ments  will  integrate  well  with  the  work  on  the  recirculating  rotation  sensor 
since  the  experimental  facility  being  built  is  ideally  suited  for  testing 
some  of  these  new  ideas.  The  second  major  problem  related  to  the  ultimate 
sensitivity  of  fiber  rotation  sensors  arises  from  the  effects  of  birefrin¬ 
gence  and  the  coupling  between  the  two  orthogonal  polarizations  in  a  single 
mode  fiber.  We  had  previously  pointed  out  the  importance  of  polarization 
and,  in  particular,  the  importance  of  reciprocity  in  operation  of  fiber 
rotation  sensors  and  we  have  now  developed  a  mathematical  formalism  which 
can  treat  birefringence  effects  in  a  most  general  way,  including  environmen¬ 
tal  effects  which  are  important  both  for  the  Rayleigh  scattering  and  for  the 
cross  talk  between  two  polarization  channels.  As  a  result  of  this  theoreti¬ 
cal  analysis  we  have  developed  the  idea  that  a  fiber  rotation  sensor  using 
unpolarized  light,  that  is,  with  no  correlation  between  the  two  orthogonal 
polarization  channels,  could  lead  to  greater  sensitivity  in  a  gyro.  In  the 
body  of  the  report  the  details  of  the  formalism  we  have  developed  to  treat 
birefringent  effects  and  the  measured  properties  of  such  an  unpolarized 
gyro  are  described.  These  ideas  are  being  pursued  further  under  the  con¬ 
tinuing  contract. 
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II..  PASSIVE  RECIRCULATING  SAGNAC  SYSTEM 


A.  SUMMARY  OF  PROGRESS  DURING  REPORTING  PERIOD 

The  work  on  the  passive  recirculating  Sagnac  (PRS)  system  during  the 
reporting  period  has  been  concerned  with  design  and  construction  of  a  new 
system.  This  new  system  features  the  use  of  all-fiber  components  in  the 
sensing  loop  and  interferometer. 

The  new  design  is  based  on  the  results  obtained  under  the  program 
during  the  preceding  year,  which  validated  the  basic  mechanisms  of  the  PRS 
system.  That  system  employed  hybrid,  bulk-optic  components  and  unguided 
optical  beams  in  space  for  all  portions  of  the  system  except  for  the  sensing 
loop  proper.  As  shown  in  the  report  for  that  period^  that  first  system,  in 
addition  to  serving  the  purpose  of  establishing  the  basic  principle,  also 
demonstrated  clearly  the  non-feasibility  of  pursuing  recirculating  systems 
further  with  bulk  optics  components. 

In  addition  to  their  obvious  size,  cost  and  stability  disadvantages, 
bulk  components  also  suffer  from  multimoding  and  alignment  problems,  and 
the  experience  with  the  first  system  made  it  clear  that  simultaneous  align¬ 
ment  for  efficient  injection,  balanced  counter-propagating  wave  amplitudes, 
and  low  loss  recirculation,  involved  mutually  interacting  adjustments  which 
were  extremely  difficult  to  program,  carry  out,  and  stabilize.  It  also  be¬ 
came  clear  that  all-fiber,  all-single  mode  components  should  completely  cir¬ 
cumvent  all  of  these  problems. 

Vinal  Report  for  Grant  No.  AFOSR-76-3070,  "Research  on  New  Approaches  to 
Optical  Systems  for  Inertial  Rotation  Sensing,"  Edward  L.  Ginzton  Labora¬ 
tory  Report  No.  3094,  Stanford  University  (February  1980). 


The  bulk-optics  version  was  constructed  at  a  time  when  no  all -fiber 
components  existed.  The  need  for  such  components  for  the  PRS  system  was 
the  initial  driving  force  for  a  program  here  which  resulted  in  the  develop¬ 
ment  of  components  of  this  kind.  It  was  realized  that  such  components  would 
have  wide  applicability  to  other  programs  both  here  and  elsewhere,  for  gen¬ 
eral  signal  processing  devices  using  optical  fibers.  The  first  of  such 
single-mode  fiber  devices,  a  directional  coupler,  was  developed  under  Joint 
Services  support  (JSEP)  and  this  was  followed  by  subsequent  development  of 
a  polarizer,  which  used  the  basic  coupler  format,  and  a  polarization  control¬ 
ler,  both  under  industrial  support  from  Atlantic  Richfield  Company.  The 
designs  and  fabrication  techniques  for  these  components  are  now  established 
and  are  available  to  the  present  program  for  use  in  an  all-fiber  PRS  system. 

During  the  reporting  period,  we  suffered  from  a  materials  problem.  We 
have  switched  from  a  jacketed  fiber  (ITT)  used  during  the  preceding  year^  to 
a  new  fiber  (Corning).  The  purpose  was  to  change  from  the  visible  light  (HeNe) 
used  in  the  bulk-optic  system  to  near  IR  light  (Nd;YAG).  This  change  had  the 
advantages  of  lower  propagation  loss  for  the  sensing  loop  and  direct  compati¬ 
bility  with  the  Nd:YAG  single  crystal  fiber  amplifier  being  developed  under 
another  program,  and  which  will  be  used  to  convert  the  PRS  system  into  an 
active  recirculating  system  (ARS). 

The  Corning  fiber,  although  having  the  desired  transmission  characteris¬ 
tics  at  1.06  M  wavelength,  was  unjacketed.  No  jacketed  fiber  designed  for 
operation  at  that  wavelength  was  available  at  that  time.  During  the  report¬ 
ing  period  we  did  fabricate  operational  all-fiber  components  using  that  fiber. 
However,  the  fiber  was  found  to  be  very  difficult  to  work  with.  The  procedures 
which  were  productive  with  jacketed  fiber  were  not  satisfactory  for  the  un¬ 
jacketed  fiber,  and  the  yield  was  low.  Also,  fiber-to-fiber  splices  using 
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this  fiber  showed  high  insertion  loss  (several  dB  per  splice),  and  the 
fabrication  problems  with  this  fiber  made  it  non-feasible  to  attempt  to 
fabricate  multiple  components  on  an  unbroken  fiber  to  avoid  splices. 
Multiple  components  were  successfully  constructed  in  this  way  using  the 
ITT  jacketed  fiber,  fo^  operation  at  .9  p.  This  approach  has  important 
advantages  for  the  PRS  system  because  of  the  sensitivity  to  splice  losses 
when  the  recirculating  waves  make  multiple  transits  across  the  splices. 

By  the  time  this  stage  was  reached,  a  jacketed  fiber  designed  for 
operation  in  the  1.06  p  range  became  commercially  available  (Corning) 
and  the  decision  was  then  made  to  change  to  the  new  fiber.  This  occurred 
near  the  end  of  the  reporting  period,  and  procurement  was  initiated  at 
that  time,  with  expectation  of  delivery  in  the  first  few  months  of  the 
following  period. 

During  the  reporting  period  the  design  for  the  all-fiber  PRS  system 
was  undertaken  and  completed.  By  the  end  of  the  period,  techniques  for 
fabrication  of  all  components  needed  for  the  fiber  portions  of  the  system 
were  on  hand,  and  a  new  fiber  was  on  order  which  was  expected  to  allow 
fabrication  of  these  components  on  a  single  unbroken  strand  of  fiber.  De¬ 
sign,  procurement  and  construction  activities  for  all  other  (non-fiber) 
portions  of  the  new  PRS  system  continued  during  the  reporting  period. 

B.  SYSTEM  DESIGN 

1 .  Basic  Recirculating  System 

In  the  basic  recirculating  Sagnac  system,  shown  in  block  form  in 
Fig.  1,  light  from  a  pulsed  source  is  sent,  via  an  optical  junction,  into 
a  multiturn  fiber  sensing  loop,  where  it  makes  multiple  circulations  around 
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PULSE 

GENERATOR 


the  loop.  The  light  recirculating  in  the  loop  is  nondestructi vely  sampled 
by  the  junction  and  sent  to  an  optical  detector  and  associated  video  signal 
processor  and  display  which  measure  inertial  rotation  of  the  sensing  loop 
about  an  axis  perpendicular  to  its  plane. 

The  junction  has  the  property  that  light  entering  at  A  is  split  into 
two  equal-amplitude  parts  at  B  and  C  ,  leading  towaves  propagating  clock¬ 
wise  and  counterclockwise,  respectively,  around  the  loop.  It  also  has  the 
property  of  closing  the  loop  upon  itself,  so  that  light,  once  propagating 
around  the  loop,  continues  to  do  so  in  multiple  transits.  The  junction  also 
has  the  property  that  it  couples  a  given  portion  of  both  counter-propagating 
light  waves,  recirculating  in  the  loop,  back  out  of  the  loop  and  sends  them 
to  the  detector. 

The  phases  of  the  light  waves  are  measured  after  each  transit  of  the 
waves  around  the  loop.  To  facilitate  this,  the  output  of  the  laser  is  pulsed 
with  a  pulse  width  somewhat  less  than  the  optical  transit  time  once  around 
the  complete  loop.  The  loop  is  active,  containing  an  optical  amplifier  and 
other  elements  (G  in  the  figure)  to  conserve  pulses  circulating  in  the  loop. 
The  amplifier  compensates  for  optical  attenuation  due  to  propagation  loss  in 
the  loop,  energy  extracted  by  the  junction  in  the  sampling  process,  and  any 
internal  loss  in  the  junction.  This  leads  to  a  time  sampled  output  waveform 
whose  envelope  is  sinusoidal  in  time.  Measurement  of  the  frequency  of  this 
waveform  provides  a  measure  of  the  angular  velocity  of  the  system.  This  is 
an  essentially  linear  system,  and  as  such  is  free  of  mode  locking  phenomena 
experienced  with  ring  laser  gyros. 

Consider  a  single  pulse  extracted  from  the  laser,  leading  to  two  pulses 
traveling  in  the  loop,  one  clockwise  and  the  other  counterclockwise.  After 
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these  pulses  have  travelled  once  around  the  N-turns  of  the  loop,  they  are 
reintroduced  into  the  loop  for  a  second  transit,  and  this  process  is  re¬ 
peated  until  a  total  of  K  transits  have  been  completed.  The  amplifier 
in  the  loop  allows  a  large  number  of  recirculations  K  ,  for  large  inte¬ 
gration  time  against  the  rotation,  and  sensitivity  to  small  rotation  rates. 

Each  time  the  pulses  complete  a  transit  around  the  loop,  they  are 
sampled,  and  these  samples  are  combined  and  sent  to  the  detector.  As  the 
original  pulses  continue  to  recirculate  around  the  loop,  they  generate  a 
train  of  pulse  samples  at  the  detector,  as  a  function  of  time.  It  is  im¬ 
portant  to  emphasize  that  this  entire  train  of  output  pulses  results  from 
one  single  pulse  extracted  from  the  laser.  The  effect  of  rotation  is  to 
shift  the  phases  of  the  optical  pulses  circulating  around  the  loop.  The 
clockwise  and  counterclockwise  pulses  are  shifted  in  phase  in  opposite  di¬ 
rections  by  the  rotation,  so  that  after  exposure  to  rotation  they  possess 
a  difference  in  phase,  or  relative  phase  shift.  This  relative  phase  shift 
is  proportional  to  the  rotation  rate  and  also  proportional  to  the  length 
of  time  the  pulses  have  been  exposed  to  the  rotation.  As  the  pulses  con¬ 
tinue  to  circulate  around  the  loop,  they  accumulate  phase  shift  Aj3(t) 
which  increases  linearly  with  time  t  .  The  detector  is  a  phase  detector 
whose  output  is  sinusoidal  in  the  argument  A0(t)  .  That  is,  for  a  fixed 
rotation  rate,  the  output  of  the  detector  is  a  sinusoidal  wave  in  time. 

The  frequency  of  this  sinusoidal  waveform  is  the  same  as  the  beat  frequency 
in  a  standard  ring  laser  gyro  (RLG). 

We  refer  to  the  system  just  described  as  an  active  recirculating  Sagnac 
system  (ARS)  because  of  the  presence  of  the  amplifier  in  the  sensing  loop. 
When  operated  as  a  passive  recirculating  system  (PRS)  without  the  amplifier, 
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the  output  pulse  train  is  a  damped  sinusoid  whose  amplitude  decreases 
experimentally  with  time. 

2 .  All -Fiber  Interferometer 

The  design  for  the  optical  junction  in  single-mode  fiber  form  is 
shown  in  Fig.  2(b),  where  the  terminals  A,  B,  C,  D  correspond  to  those 
in  Fig.  1.  DCl  and  DC2  are  single-mode  fiber  directional  couplers. 

DCl  is  a  3  dB  coupler  which  splits  the  light  pulses  entering  at  A  into 
two  equal  amplitude  pulses  at  E  and  F,  which  enter  coupler  DC2  at  oppo¬ 
site  ends.  DC2  is  a  directional  coupler  with  a  coupling  coefficient 
(coupling  between  B  and  C  and  between  E  and  F)  near  unity.  Thus  a  large 
fraction  of  the  light  at  E  and  F  returns  to  F  and  E,  respectively,  while 
a  small  fraction  enters  the  loop  to  form  counterclockwise  and  clockwise 
traveling  pulses,  respectively.  Pulses  once  in  the  loop  remain  in  the 
loop  by  crossing  from  one  side  of  DC2  to  the  other  following  the  path 
shown  by  the  dashed  arrow,  with  low  insertion  loss  because  of  the  high 
coupling  coefficient  between  the  two  halves  1  and  2  of  the  coupler,  and 
the  low  internal  dissipative  loss  of  the  coupler. 

Making  high  quality  splices  on  single  mode  fibers  is  difficult  be¬ 
cause  of  the  small  core  diameter.  Commercially  available  multimode  spli¬ 
cing  equipment  yields  splice  losses  of  1-3  dB  when  used  on  single  mode 
fibers.  Single  mode  splicing  equipment  with  splice  losses  of  around  0.1 
dB  is  not  currently  available.  In  order  to  determine  base-line  perfor¬ 
mance  for  recirculating  systems,  the  junction  of  Fig.  2  is  designed  to 
be  fabricated,  together  with  the  sensing  loop,  with  no  fiber  splices  re¬ 
quired.  This  calls  for  coupler  halves  1,  2,  3,  and  4  to  be  fabricated 
directly  on  the  ends  of  the  single-mode  fiber  sensing  loop  as  at  (a)  and 
then  assembled  as  at  (b). 
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1,  2,  3,  4:  DIRECTIONAL  COUPLER  HALVES 


1  4 


FIBER 

SENSING 

LOOP 


2  3 

(a)  UNASSEMBLED 

JUNCTION 


(b)  ASSEMBLED 


Fig.  2--Schematic  for  spliceless  sensing 
loop  and  junction  assembly. 
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3.  Overall  System 


A  diagram  of  the  overall  PRS  system  is  given  in  Fig.  3.  The  system 
shown  is  passive,  without  an  amplifier  in  the  sensing  loop,  pending  com¬ 
pletion  of  a  Nd:YAG  fiber  amplifier  which  is  now  under  development  in  this 
laboratory  under  other  auspices  and  will  be  added  at  a  later  date. 

The  optical  source  consists  of  a  Nd;YAG  laser  operating  at  1.06  u. 
telescope,  acousto-optic  beam  deflecting  modulator,  and  beam  steering  device 
leading  up  to  the  fiber  input. 

Polarization  controllers  PCI  and  PC2  are  necessary  to  compensate  the 
system  for  the  birefringence  in  the  fiber,  in  order  to  balance  the  inter¬ 
ferometer  at  zero  rotation  rate. 

Adjustments  of  the  fiber  circuit  are  expected  to  be  simple  and  straight¬ 
forward  in  sharp  contrast  to  earlier  models  using  bulk  optic  components  with¬ 
in  the  interferometer.  There  are  no  optical  alignments  to  be  made  within 
the  interferometer. 


C.  COMPONENT  DEVELOPMENT 


The  status  of  the  various  components 
cussed  in  this  subsection. 

1 .  Interferomter  Components 
Sensing  Loop.-  The  sensing  loop  will 
15  cm  and  will  be  wound  with  fiber  having 

Fiber  type:  single 


in  the  circuit  of  Fig.  3  is  dis¬ 


have  a  diameter  of  the  order  of 
the  following  specifications: 
mode 


Length  =  1  km 
Core  diameter  =  6  p 
Cladding  diameter  =  125  y 
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DEVICE 
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-Schematic  for  PRS  system  with 
all -fiber  interferometer. 


Optical  attenuation  (1.06  u)  =  2  dB/km 
Jacket  diameter  =250  y  . 

The  loss  in  modern  optical  fibers  is  limited  by  Rayleigh  scattering 
which  varies  inversely  with  the  fourth  power  of  the  wavelength.  To  de¬ 
crease  this  loss  requires  use  of  wavelengths  in  the  near  infrared.  A 
wavelength  of  1.06  ym  was  chosen  because  of  the  availability  of  powerful 
sources,  the  availability  of  fiber  with  loss  of  around  2  dB/km,  and  com¬ 
patibility  with  a  Nd:YAG  amplifier  under  development.  The  sensing  loop 
length  was  chosen  to  be  1  km  because  of  the  availability  of  fiber  in  that 
length  and  the  high  rotation  sensitivity  which  it  provides,  and  because 
with  this  length  of  fiber,  pulses  of  several  ys  duration  can  be  used  which 
simplifies  the  equipment  requirements  for  generation  and  detection  of  pulses. 

Directional  Couplers.  A  single-mode  fiber-to-fiber,  evanescent  field 

2 

directional  coupler  has  been  developed  in  this  laboratory,  as  noted  above, 
and  this  design  will  be  followed  in  fabricating  the  couplers  DCl  and  DC2 
on  the  fiber  just  described.  Briefly,  in  this  coupler  design,  a  fiber  is 
bonded  into  a  slot  cut  in  a  quartz  substrate  (Fig.  4a).  The  face  of  the 
block  is  polished  to  remove  a  significant  portion  of  the  cladding,  thus 
allowing  access  to  the  evanescent  fields  of  the  guided  mode  (Fig.  4b). 

When  two  such  substrates  are  placed  in  contact  (Fig.  4c),  strong  coupling 
exists  between  the  two  fibers.  This  coupling  is  highly  directional  and 
can  result  in  complete  power  transfer  in  approximately  1  millimeter  of 
distance.  The  coupling  can  be  varied  by  transverse  displacement  of  the 
fibers  (Fig.  4d). 

^R.  A.  Bergh,  G.  Kotler,  and  H.  J.  Shaw,  "Single-Mode  Fibre  Optic  Direc¬ 
tional  Couplers,"  Electronics  Letters  7,  260  (27  March  1980). 
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CLADDING 


Figl  4--Fiber  to  fiber,  evanescent  field  coupler,  (a)  Fiber  bonded 
into  slot  cut  in  fused  silica  substrate,  (b)  Cladding  par¬ 
tially  removed  to  expose  evanescent  field,  (c)  Assembled 
coupler  side  view,  (d)  Tuning  coupler  (end  view)  by  trans¬ 
verse  displacement  of  blocks. 
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This  coupler  design  is  characterized  by  mechanically  variable  coup¬ 
ling  coefficient,  high  directivity,  low  internal  loss  and  low  backscat- 
tering.  Coupler  DC-1  will  be  constructed  with  maximum  coupling  coeffi¬ 
cient  of  70%,  which  can  then  be  set  experimentally  to  50%  (3  dB  coupler) 
during  operation.  Coupler  DC2  will  be  constructed  with  slight  over 
coupling  at  the  maximum  setting  (interaction  length  >  1/2  wavelength), 
which  can  then  be  set  experimentally  to  any  desired  value  during  opera¬ 
tion.  In  the  PRS  system,  there  is  an  optimum  value  for  this  coupling 
coefficient  k  which  maximizes  the  number  n  of  observable  pulses  in 
the  output  pulse  train,  as  n  is  given  by 


n 


£n  —  +  £n 


Pg  (1  -  k)‘ 


£N(kR) 


(1) 


where 

Pq  =  input  power  level 

P  =  minimum  detectable  power  level 
m 

R  =  power  transmission  coefficient  of  sensing  loop  = 
exp(-a£)  where  a  is  the  power  attenuation  coef¬ 
ficient  of  the  fiber  and  £  is  the  fiber  length. 

For  example,  the  optimum  value  of  k  varies  from  .975  to  .75  as  the 
fiber  loss  varies  from  0  to  3  dB,  and  the  corresponding  value  of  n 
varies  from  15  to  7.  Thus  a  coupler  which  is  capable  of  slight  over¬ 
coupling  can  be  tuned  through  this  entire  range,  and  set  experimentally 
to  match  the  attenuation  constant  of  the  particular  fiber  used,  so  as 
to  maximize  the  number  of  recirculations. 
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For  coupler  DC2  the  low  internal  loss  in  the  present  coupler  design 
(measured  <  0.5  dB)  is  crucial,  as  it  is  directly  in  series  with  the  pulse 
recirculation  path  in  the  loop. 

Polccpization  Controllers.  The  polarization  controllers  PCI  and  PC2 
(Fig.  3)  are  devices  which  can  transform  any  state  of  polarization  (in 
general  elliptic  with  arbitrary  ellipticity  and  axis  orientation)  at  their 
input  end  to  any  other  state  of  polarization  at  the  output  end.  A  design 

3 

for  this  device  has  been  developed  under  another  program  in  this  labora¬ 
tory.  These  controllers  do  not  require  any  fiber  splices,  and  do  not  in¬ 
volve  any  surgery  on  either  the  fiber  or  its  jacket.  They  involve  two 
small -diameter  (1  or  2  cm  depending  on  fiber  characteristics)  fiber  loops 
whose  bend-induced  birefringence  constitutes  the  equivalent  of  fractional- 
wave  plates  of  very  low  loss,  whose  principle  axes  can  be  rotated  experi¬ 
mentally  to  produce  any  desired  transformation  of  the  state  of  polarization. 
2.  External  Components 

This  section  discusses  the  components  in  Fig.  3  outside  of  the  fiber 
circuit. 

Laser.  The  external  source  is  a  General  Photonics  Model  TWO-15,  Ne:YAG 

laser  operating  at  1.06  ym.  The  power  output  is  2  watts  in  the  TEM^^  mode 

2 

with  a  polarization  of  1000:1.  The  1/e  beam  diameter  is  21  mm.  The  laser 
is  equipped  with  a  closed  loop  water-to-air  heat  exchanger,  thus  permitting 
the  mounting  of  the  entire  laser  system  on  the  rotation  table.  The  amplitude 
stability  of  the  laser  was  measured  at  around  10%. 

Pulse  Modulator.  An  important  design  consideration  is  the  requirement 
for  a  large  number  of  recirculations.  The  ultimate  limit  on  recirculations 

3 

H.  C.  Lefevre,  "Single-Mode  Fibre  Fractional  Wave  Devices  and  Polarisation 
Controllers,"  Electronics  Letters  16,  20,  778  (25  September  1980). 
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is  due  to  noise.  In  a  poorly  designed  system,  a  significantly  lower 
limit  on  recirculations  can  result  from  sensing  loop  round  trip  loss 
and  finite  modulation  depth  of  the  optical  pulse  generator.  The  non¬ 
zero  power  transmitted  by  the  pulse  generator  in  the  interpulse  period 
results  in  an  optical  background  into  which  the  recirculations  decay 
because  of  sensing  loop  loss.  The  expected  number  of  recirculations, 
n  ,  is  approximately  given  as 

n  =  d/ail  (2) 

where  d  is  the  modulation  depth  and  a£  is  the  total  round  trip  loss 
in  the  sensing  loop.  For  large  n  ,  Eq.  (2)  requires  large  modulation 
depths  and/or  low  round  trip  losses.  The  round  trip  loss  in  the  sensing 
loop  has  four  components:  loss  in  the  optical  fiber,  splice  losses,  in¬ 
sertion  loss  of  coupler  DC2,  and  sampling  loss.  Sampling  loss  is  due  to 
coupling  a  nonzero  amount  of  optical  power  from  the  sensing  loop  each 
recirculation  for  measurement  of  the  accumulated  Sagnac  phase  shift. 

Based  on  the  above  considerations,  the  choice  of  an  optical  modula¬ 
tor  was  largely  concerned  with  modulation  depth  capability.  Q-switching 
the  laser  gives  large  modulation  depths  but  also  results  in  high  peak 
powers  which  are  not  appropriate  for  use  with  single  mode  fiber  operating 
in  the  linear  range.  In  addition,  the  Q-switched  pulses  are  narrow,  which 
would  require  sophisticated  detection  equipment.  An  external  modulator  on 
a  cw  laser  can  give  long  pulses  with  low  peak  powers.  Standard  commercial 
modulators  have  modulation  depths  of  20-30  dB.  For  these  modulators,  Eq. 
(2)  predicts  6  to  10  observable  recirculations.  Custom  built  systems 
could  achieve  around  45  dB  of  modulation  depth  but  at  a  great  increase  in 
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price.  Furthermore,  this  modulation  depth  is  not  as  high  as  desired  be¬ 
cause  it  would  limit  the  number  of  recirculations  to  15  or  less. 

The  pulse  modulator  shown  in  Fig.  5  was  developed  to  provide  modu¬ 
lation  depths  of  better  than  90  dB.  The  major  components  are  an  acousto- 
optic  beam  deflector,  a  coupling  lens,  and  the  optical  fiber  which  func¬ 
tions  as  an  apodized  pinhole.  This  modulator  functions  as  follows.  When 
the  acousto-optic  beam  deflector  is  on,  the  optical  beam  is  deflected 
through  an  angle  (9)  ,  and  the  deflected  beam  is  focused  into  the  opti¬ 
cal  fiber  by  the  coupling  lens.  The  undeflected  beam  is  focused  a  distance 
A  away  from  the  core  of  the  fiber  given  by 

A  =  f9  (3) 

where  f  is  the  focal  length  of  the  coupling  lens.  For  the  values  shown 
in  Fig.  5,  A  is  106  pm.  The  spot  size  of  the  focused  beam  is  around  3 
pm,  hence  very  little  power  from  the  undeflected  beam  is  coupled  into  the 
fiber.  The  limit  on  the  modulation  depth  is  due  to  reflections  from  the 
surfaces  of  the  coupling  lens,  and  is  95  dB  for  an  uncoated  lens.  For  an 
AR  coated  lens  the  modulation  depth  can  be  increased  to  about  115  dB. 
Measurements  of  the  modulation  depth  placed  it  at  better  than  46  dB.  This 
measurement  was  limited  by  the  sensitivity  of  the  measurement  system. 

Redueing  Telescope.  A  reducing  telescope  of  reduction  ratio  2.5  is 
used  to  match  the  laser  beam  to  the  acoustic  aperture  of  the  acousto-optic 
beam  deflector. 

Beam  Deflector.  A  Coherent  Model  304  acousto-optic  beam  deflector 
is  used.  The  deflection  angle  is  7.1  mrads  and  the  deflection  efficiency 
was  measured  at  55%  at  a  wavelength  of  1.06  pm.  The  transmittance  of  the 
reducing  telescope  and  beam  deflector  combination  was  measured  to  be  70% 
at  1.06  pm. 
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.  5--Optical  pulse  generator  with  large  modulation  depth 


steering  Device.  A  model  605  beam  steering  device  by  Newport  Re¬ 
search  Corporation  is  used  to  steer  the  laser  beam  onto  the  coupling 
lens.  The  transmittance  of  this  device  is  99%  at  1.06  pm. 

Coupling  Lens.  A  plano-convex  lens  of  15  mm  focal  Inegth  is  used. 

The  coupling  efficiency  into  the  fiber  was  measured  at  25%. 

Photodiodes,  Four  types  of  photodiodes  are  available:  A  Centronics 
model  OSD-1 -4  IR  enhanced  silicon  diode,  a  Centronics  model  APD-05-4,  IR 
enhanced  silicon  avalanche  diode,  a  Judsen  Infrared  model  0-16-5  germanium 
diode,  and  a  Ford  Aerospace  model  L4520  low  noise  germanium  diode. 

Transimpedance  Amplifiers.  Three  models  are  available  from  Judson 
Infrared:  model  070,  gain  50  dB,  bandwidth  5  Hz-10  MHz;  model  470,  gain 
26  dB,  bandwidth  200  Hz-30  MHz,  model  490,  gain  26  dB,  bandwidth  200  Hz- 
100  MHz. 

Pulse  Generator.  An  HP  model  801 3B  pulse  generator  is  used  for  the 
signal  input  into  the  acousto-optic  driver.  The  pulse  height  is  5V,  the 
pulse  width  2  ys  and  the  pulse  period  200  ys.  The  pulse  generator  also 
provides  a  trigger  signal  for  the  oscilloscope. 

Oscilloscope.  An  HP  model  1742A,  100  MHz  oscilloscope  is  used  to 
display  the  measured  pulse  heights  versus  time.  The  bandwidth  of  the 
oscilloscope  is  user  limited  to  20  MHz  to  reduce  noise. 
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III.  NEW  TECHNIQUES  FOR  FIBER  GYROS 


A.  INTRODUCTION 

All  existing  fiber  gyros  are  based  on  the  same  mechanism,  the  Sagnac 
effect,  in  which  EM  waves  traversing  a  closed  path  experience  a  phase 
shift  proportional  to  the  rate  of  rotation  of  the  path  in  inertial  space. 

All  such  gyros  also  use  a  multiturn  sensing  loop  for  the  closed  path,  to 
multiply  the  phase  shift  and  sensitivity.  The  various  types  of  gyros  differ 
only  in  the  details  of  processing  the  optical  signals  inside  and  outside  of 
the  sensing  loop,  and  the  electrical  signals  from  the  detector,  to  emphasize 
various  characteristics,  such  as  sensitivity,  dynamic  range,  response  time, 
linearity,  etc. 

The  sensitivity  of  fiber  gyros  is  limited  by  instabilities  and  noise, 

of  which  quantum  noise  in  the  detector  represents  the  irreducible  minimum. 

As  of  the  beginning  of  the  reporting  period,  noise  from  other  sources  in 

excess  of  this  minimum  set  the  sensitivity  limit  of  all  existing  fiber  gyros. 

We  consider  first  effects  of  polarization.  Earlier  work  under  the  pre- 
4  5 

sent  program  and  elsewhere  had  shown  that  reciprocity  was  a  basic  require¬ 
ment  for  stable  operation  of  fiber  gyros  of  all  types,  i.e.,  that  input  and 
output  polarizations  of  both  counterpropagating  pulses  should  be  the  same. 
This  was  the  most  significant  advance  since  the  inception  of  the  fiber  gyro, 
but  it  was  clear  that  the  subject  was  not  closed,  because  only  the  case  of 
linear  polarization  had  been  dealt  with.  Thus  we  turned  to  the  development 

4 

H.  Arditty,  H.  J.  Shaw,  M.  Chodorow,  and  R.  Kompfner,  Proc.  Soc.  Photo-Opt. 
Instrum.  Eng.  1^,  138  (1978). 

^R.  Ulrich  and  M.  Johnson,  Opt.  Lett.  4,  5,  152  (1979). 


-  24  - 


of  a  formalism  which  would  allow  treating  the  general  case  of  arbitrary 
polarization,  as  described  in  Sections  III.B  and  C. 

Much  of  the  problem  of  increasing  the  sensitivity  of  fiber  gyros  is 
concerned  with  reducing  noise  resulting  from  residual  environmental  ef¬ 
fects  such  as  temperature  changes  and  vibration  which  remain  when  the 
gyro  is  configured  to  be  optically  reciprocal.  It  is  important  to  empha¬ 
size  that  every  so-called  single  mode  fiber  is  birefringent  and  actually 
supports  two  polarization  modes.  The  birefringent  effects  are  dependent 
upon  environmental  effects,  and  this  leads  to  some  of  the  instabilities 
and  noise  observed  in  fiber  gyros.  Sections  III.B  and  C  deal  with  the 
problem  of  handling  these  effects  analytically.  This  analysis  forms  a 
very  important  guide  to  the  process  of  reducing  the  noise  and  increasing 
the  sensitivity  of  fiber  gyros  in  general.  This  also  has  resulted  in  a 
proposal  which  is  described  below  for  a  new  form  of  gyro,  using  unpolarized 
light  to  achieve  certain  performance  advantages. 

Second,  we  consider  the  effects  of  light  scattering.  Initial  theore¬ 
tical  work  carried  out  under  the  present  contract  showed  that  coherent 
backscattering  from  random  inhomogenieties  in  the  fiber  (Rayleigh  scattering) 
is  capable  of  reducing  the  sensitivity  of  fiber  gyros  to  a  level  several 
orders  of  magnitude  less  than  that  imposed  by  photon  noise,  contrary  to 
prior  theoretical  predictions.  A  discussion  of  this  problem,  including 
numerical  estimates  of  the  magnitude  of  the  effect  and  general  approaches 
for  reducing  it,  is  included  below  in  Section  III.D. 

It  is  important  to  emphasize  that  Rayleigh  scattering,  as  referred 
to  above,  is  a  property  of  all  fibers  in  fixed  amount,  depending  on  the 
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wavelength.  In  Section  III.D  it  is  shown  that  even  though  this  scattering 
is  coherent,  it  leads  to  direct  phase  errors  which,  through  the  effects 
dealt  with  in  Section  III.B,  are  environmentally  dependent,  thus  are  non- 
systematic  and  cannot  be  calibrated,  and  must,  therefore,  be  reduced  in 
magnitude  to  acceptable  levels.  This  work  is  expected  to  have  a  major 
impact  in  increasing  the  sensitivity  of  fiber  gyros  in  general. 

The  analyses  of  Sections  III.B,  C  and  D  are  applied  to  single  pass 
gyros  using  cw  light.  This  is  the  simplest  case  and  was  the  appropriate 
starting  point.  It  has  led  to  results  of  direct  importance  to  practical 
gyros  of  this  type  which  were  under  development  in  many  laboratories. 
Extensions  of  this  work  to  recirculating  gyros  of  the  type  described  in 
Section  II  will  be  treated  in  later  reports. 

B.  BIREFRINGENCE  EFFECTS 

1 .  Introduction 

The  single  mode  optical  fibers  used  in  fiber  optic  gyroscopes  are 
birefringent  because  they  propagate  two  modes  which  have  slightly  different 
propagation  constants.  These  two  guided  modes  are  the  orthogonal  polariza¬ 
tions  of  the  HEii  mode.  In  an  ideal  fiber,  the  propagation  constants  of 
these  modes  would  be  identical.  Small  asymmetries  enter  into  the  fiber 
during  the  manufacturing  process  causing  the  fiber  to  be  birefringent. 

Intimately  associated  with  birefringence  is  the  polarization  of  the 
light  which  propagates  through  a  birefringent  medium.  To  specify  the 
polarization,  the  state  and  degree  or  polarization  must  be  specified. 

The  effect  of  propagation  through  a  birefringent  medium  is  to  change  the 
state  of  polarization  but  not  the  degree  of  polarization. 
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The  general  state  of  polarization  is  elliptic,  and  there  are  two 
familiar  limiting  cases,  namely  linear  and  circular.  These  states  are 
shown  in  Fig.  6.  The  figures  are  produced  by  the  tip  of  the  electric 
field  vector  as  it  moves  in  a  plane  perpendicular  to  the  propagation 
direction.  A  complete  description  of  the  state  of  polarization  includes 
ellipticity,  the  orientation  with  respect  to  some  coordinate  system,  and 
the  direction  (clockwise  or  counterclockwise)  in  which  the  electric  field 
vector  rotates. 

The  degree  of  polarization  is  a  statistical  description  of  how  the 
state  of  polarization  evolves  in  time.  If  the  degree  of  polarization  is 
unity,  then  the  state  of  polarization  is  constant  in  time.  If  the  degree 
of  polarization  is  zero,  then  the  state  of  polarization  changes  randomly 
in  time  through  all  possible  states. 

Single  pass,  fiber  optic  gyroscopes  measure  rotation  rate  by  measuring 
the  optical  intensity  of  the  recombined  counter-propagating  light  waves. 

This  intensity  is  a  function  of  the  phase  shift  between  the  counter-propa¬ 
gating  light  waves.  Ideally,  this  phase  shift  equals  the  Sagnac  phase  shift 
which  is  proportional  to  the  rotation  rate.  The  birefringence  of  the  opti¬ 
cal  fiber  can  introduce  additional  phase  shifts  which  may  be  indistinguish¬ 
able  from  the  rotation  induced  phase  shift.  These  birefringence  dependent 
phase  shifts  are  not  constant  because  of  the  extreme  environmental  sensiti¬ 
vity  of  the  fiber  birefringence. 

If  polarization  selective  components  exist  in  the  fiber  gyroscope, 
then  a  change  in  the  optical  intensity  results  from  a  change  in  the  polari¬ 
zation  state  due  to  fiber  birefringence.  This  birefringence  dependent  in¬ 
tensity  change  may  be  indistinguishable  from  a  rotation  induced  intensity 
change. 
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Fig.  6--Electric  field  trajectories  in  a  plane 
normal  to  propagation  direction  for 
various  states  of  polarization. 
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The  net  result  of  these  birefringence  dependent  phase  and  intensity 
changes  can  be  large  errors  in  the  measurement  of  rotation  rates.  These 
birefringence  dependent  changes  must  be  eliminated  or  reduced  in  practical 
fiber  gyroscopes,  and  this  requires  an  understanding  of  the  interaction  of 
polarization  and  fiber  birefringence  in  these  devices. 

2.  Theoretical  Development 

A  mathematical  formalism  has  been  developed  which  allows  calculation 
of  the  response  of  a  single  pass,  fiber  gyroscope  for  any  input  polariza¬ 
tion  and  any  fiber  birefringence.  This  formalism  is  both  a  generalization 

5 

and  refinement  of  an  approach  reported  by  Ulrich  and  Johnson  which  was 
limited  to  input  polarizations  whose  degree  of  polarization  is  unity,  i.e., 
where  the  state  of  polarization  is  constant  in  time. 

The  state  of  polarization  will  be  described  in  terms  of  basic  states. 
The  basic  states  are  chosen  as  the  x  and  y  linear  polarization  states. 
Any  polarization  state  can  be  expressed  as  a  weighted  sum  of  these  base 
states.  The  degree  of  polarization  is  described  by  the  correlation  of  the 
coefficients  of  these  base  states. 

The  coherence  matrix  formalism  is  chosen  to  represent  the  degree  and 
state  of  polarization.  The  coherence  matrix  J  is  defined  in  terms  of 
the  coefficients  of  the  base  states  as 

(<a^(t)  a*(t)>  <a^(t)  a*(t)>\ 

I  (4) 

<ay(t)  a*{t)>  <;ay(t) 

where  a  (t)  and  a  (t)  are  the  instantaneous  values  of  coefficients  of 
X  y' 

the  X  and  y  linear  polarization  base  states,  the  asterisk  denotes  the 
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complex  conjugate  and  the  brackets  denote  the  infinite  time  average.  Fur¬ 
ther  details  of  the  coherence  matrix  and  its  properties  are  given  in  Appen¬ 
dix  B. 

The  optical  fiber  is  modeled  as  a  general,  four  port  linear  system. 

The  model  for  this  is  shown  in  Fig.  1  of  Appendix  A.  The  four  ports  are 
the  X  and  y  linear  polarization  base  states  at  either  end  of  the  fiber. 
Also  seen  in  the  figure  are  the  four  channels  which  connect  the  four  ports 
as  shown.  The  channel  is  an  abstract  pathway  along  which  light  propagates. 
T,his  model  of  the  fiber  can  be  described  by  two  Jones  matrices,  one  for 
propagation  in  each  direction.  Ulrich  gave  the  general  form  for  these 
Jones  matrices  T  in  the  presence  of  rotation,  which  are  given  as  Eqs.  2 
and  3  of  Appendix  A. 

Explicit  forms  for  the  elements  t(ki)  have  been  derived.  This  deri¬ 
vation  is  treated  in  Appendix  C.  The  Jones  matrices  describing  the  fiber 
can  be  rewritten  as 


21  ^  ® 


'l2  ®  ^ 


cosOe’^'^^^ 

sinOe"^'^'^^ 

-sinGe"^^^^ 

cosGe 

cosGe^^'^^ 

-sinGe-^^^^ 

sinGe"'^^^^ 

cosGe""^^^^ 

(5a) 


(5b) 


In  these  matrices,  0^  is  the  Sagnac  phase  shift,  0  is  the  average  phase 
shift  along  the  fiber,  0  is  the  state-of-polarization  rotation  parameter 
(i.e.,  the  state  of  polarization  is  rotated  by  9),  C  is  the  linear  bire¬ 
fringence  of  the  fiber,  and  0  is  a  reciprocal  birefringence  dependent  phase 

^Note  that  in  Appendix  C,  notation  is  ^12  ’  ^21  ' 
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shift,  which  can  be  indistinguishable  from  the  Sagnac  phase  shift  in  certain 
cases.  The  parameters  0  ,  9  ,  C  .  and  0  depend  on  the  birefringence  in 
the  fiber  and  are  environmentally  sensitive.  These  parameters  are  defined 
more  completely  in  Appendix  C. 

A  fiber  gyroscope  typically  contains  many  components.  These  components 
include  beam  splitting  and  combining  elements,  polarizers,  polarization  con¬ 
trollers,  phase  shifters  and  phase  modulators.  Each  component  is  represented 
by  one  or  more  Jones  matrices.  The  Jones  matrix  of  a  fiber  gyroscope  can  be 
written  in  terms  of  the  Jones  matrices  of  the  components,  using  matrix  multi¬ 
plication  and  addition.  Let  M  denote  the  Jones  matrix  of  the  fiber  gyro¬ 
scope.  The  coherence  matrix  of  the  output  light  J^^^^  is  related  to  the 
coherence  matrix  of  the  input  as 


"^out 


(6) 


where  m’*'  is  the  hermitian  adjoint  of  M  .  The  output  intensity  can  be 
computed  from  the  output  coherence  matrix. 

3.  Unpolarized  Operation 

This  mathematical  formalism  was  applied  to  a  single  pass,  cw  fiber 
optic  gyroscope  operated  with  unpolarized  input  light.  The  coherence  ma¬ 
trix  for  unpolarized  light  of  unit  intensity  is 


/1/2  0 

'^unpolarized  "  I 

\  0  1/2 

This  says  that  unpolarized  light  has  equal  average  power  in  the  two  ortho¬ 
gonal  base  states  and  that  the  coefficients  of  the  base  states  are  uncor¬ 
related,  which  means  that  the  electric  fields  in  the  two  orthogonal  base 
states  are  incoherent  and  that  they  cannot  form  an  interference  pattern. 
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The  calculated  output  intensity  for  unit  input  intensity  for  the 
unpolarized  gyroscope  in  Fig.  7  is 

I  =  1/2  +  1/2  (cos^9  -  sin^Q  cos0}  cos0g  (8) 

where  0  describes  the  rotation  of  the  state  of  polarization  by  the 
fiber,  and  0  is  the  reciprocal  birefringent  dependent  phase  shift. 

These  quantities  are  derived  in  Appendix  C.  In  obtaining  Eq.  (8),  the 
beam  splitting/combining  elements  are  assumed  ideal.  This  means  that 
they  are  lossless,  polarization  insensitive,  and  have  a  50-50  splitting 
ratio. 

A  plot  of  I  versus  0^  (Sagnac  phase  shift)  for  Eq.  (8)  is  shown 
in  Fig.  8.  The  curves  in  Fig.  8  are  described  by  different  values  of  0 
and  0  which  are  functions  of  the  fiber  birefringence.  Three  features 
of  these  curves  are  of  interest.  First,  the  maxima  and  minima  of  the 
curves  always  occur  at  the  same  value  of  0^  .  No  translation  of  the 
curves  along  the  0^  axis  exists  for  changing  values  of  0  and  0  . 

If  this  type  of  translation  was  present,  then  the  Sagnac  phase  shift 
would  be  indistinguishable  from  0  ,  the  reciprocal  birefringence  depen¬ 
dent  phase  shifts.  Second,  the  average  value  of  the  intensity  for  these 
curves  is  constant.  The  maxima  and  minima  of  the  curves  expand  or  con¬ 
tract  about  the  average  value  with  varying  birefringence.  In  conventional, 
polarized,  single  pass  fiber  optic  gyroscopes  the  value  of  the  minima  is 
fixed,  while  the  value  of  the  maxima  of  the  curves  changes  with  varying 
birefringence.  Third,  all  the  response  curves  in  Fig.  8  intercept  each 
other  at  a  family  of  points  which  are  characterized  by  the  Sagnac  phase 
shift  being  an  odd  multiple  of  7r/2  radians.  These  points  of  intersection 
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F1g.  7--Schematic  of  a  basic,  reciprocal,  unpolarized  fiber  gyroscope. 
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-Normalized  output  intensity  versus  Sagnac  phase  shift  for  reciprocal, 
unpolarized  gyroscope.  The  curves  are  described  by  the  birefringent 
parameters  (6,0).  (a)  (45,0),  (b)  (20,25),  (c)  (30,15),  (d)  (40, 


of  the  response  curves  are  called  stable  operating  points.  If  the  gyro¬ 
scope  is  operating  at  one  of  these  points,  the  detected  intensity  is  in¬ 
dependent  of  changing  birefringence.  The  gyroscope  is  stable  with  respect 
to  reciprocal  changes  in  the  fiber  birefringence  at  these  points.  An  un¬ 
polarized  gyroscope  has  been  built  and  operated  during  the  reporting  period. 
It  confirmed  the  theoretical  predictions  above. 

Further  calculated  cases,  and  a  description  of  the  experimental  results, 
are  given  in  Appendix  C. 

Advantages  of  unpolarized  fiber  optic  gyroscopes  over  conventional  po¬ 
larized  fiber  optic  gyroscopes  are  expected  to  be  a  reduced  part  count  (no 
polarizer  is  required),  reduced  environmental  sensitivity,  and  a  reduction 
in  coherent  backscattered  noise  due  to  the  use  of  two  spatially  incoherent 
modes  (i.e.,  the  two  orthogonal  polarization  modes). 

C.  PROPOSED  UNPOLARIZED  GYROSCOPE 

The  design  shown  in  Fig.  9  is  proposed  for  an  unpolarized  gyroscope. 

It  takes  full  advantage  of  the  operating  characteristics  described  above 
for  unpolarized  gyroscopes,  to  provide  a  stable  output  which  is  linear  in 
the  rotation  rate  over  a  broad  range. 

The  gyroscope  is  designed  so  that  the  operating  point  (i.e.,  values 
of  output  intensity  and  phase  shift  between  the  counterpropagating  waves) 
remains  at  or  near  one  of  the  stable  operating  points.  This  is  accomplished 
by  sensing  the  difference  between  the  output  intensity  and  the  average  in¬ 
tensity,  and  using  feedback  applied  to  the  variable  nonreciprocal  phase 
shifter  in  the  sensing  loop  to  cause  the  output  intensity  to  equal  the 
average  intensity.  The  amount  of  feedback  applied  is  a  measure  of  the  ro¬ 
tation,  and  is  linear  in  the  rotation  rate. 
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Fig.  9 — Schematic  of  proposed  unpolarized  fiber  gyroscope. 
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At  rest  a  tt/2  radian  phase  shift  exists  between  the  counterpropaga- 
ting  waves  due  to  the  nonreciprocal  phase  shifter.  This  nonreciprocal 
phase  shift  0^^  is  additive  to  the  Sagnac  phase  shift  (i.e.,  0^  ->  0^  + 
0nr)-  This  biases  the  gyroscope  to  the  first  stable  operating  point.  Eq. 
(8)  becomes 

I  =  1/2  +  1/2  (cos^9  -  sin^9  cos0)  cos  (11/2)  =  1/2  .  (9) 

The  output  intensity  at  rest  equals  the  average  intensity.  When  the 
gyroscope  is  rotated,  the  Sagnac  phase  shift  is  no  longer  zero.  Before 
feedback  is  applied  to  the  nonreciprocal  phase  shifter,  Eq.  (9)  becomes 

I  =  1/2  +  1/2  (cos^9  -  sin^0  cos0)  cos(0^  +  r/2)  .  (10) 

The  output  intensity  now  differs  from  the  average  intensity.  Feedback 
is  applied  to  the  nonreciprocal  phase  shifter  so  that  the  output  intensity 
above  equals  the  average  intensity.  Eq.  (10)  becomes 

I  =  1/2  +  1/2  (cos^0  -  sin^O  cos0)  cos(0^  +  tt/2  +  0^^^)  =  1/2  •  (11) 

The  nonreciprocal  phase  shift  which  accomplishes  this  is  equal  to,  but  op¬ 
posite  in  sign  from,  the  Sagnac  phase  shift.  If  the  nonreciprocal  phase 
shifter  is  a  Faraday  cell,  then  the  drive  current  is  proportional  to  the 
nonreciprocal  phase  shift  and  hence  the  Sagnac  phase  at  the  stable  opera¬ 
ting  point.  By  measuring  the  drive  current  the  rotation  rate  is  obtained 
and  the  drive  current  is  also  linear  in  the  rotation  rate. 

The  unpolarized  gyroscope  depicted  in  Fig.  9,  implements  the  above 
scheme.  Light  from  an  unpolarized  source  is  coupled  into  the  fiber  by  a 
coupling  lens.  It  passes  through  a  reciprocity  coupler  which  is  a  fiber- 
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to-fiber  evanescent  field  directional  coupler.  Approximately  half  of  the 
light  is  sent  to  the  input  detector.  The  resulting  signal  is  used  to  com¬ 
pute  the  average  output  intensity.  This  is  done  to  compensate  for  laser 
power  and  input  coupling  fluctuations. 

The  remaining  light  propagates  to  the  splitter/combiner  coupler  where 
it  is  split  into  two  waves  of  approximately  equal  power,  which  then  propa¬ 
gate  in  opposite  directions  around  the  N-turn,  optical  fiber  sensing  loop. 
The  variable  nonreciprocal  phase  shifter  introduces  a  phase  shift  between 
the  counterpropagating  waves  in  addition  to  any  rotation  induced  phase 
shift.  A  polarization  controller  is  provided  to  reduce  any  errors  intro¬ 
duced  by  a  residual  degree  of  polarization  of  the  light  source.  It  also 
exerts  control  over  the  slope  (i.e.,  output  intensity 

at  the  stable  operating  point. 

The  light  waves  are  recombined  by  the  splitter/combiner  coupler  and 
propagate  to  the  reciprocity  coupler,  where  half  of  the  power  is  then  sent 
to  the  output  detector.  The  electrical  signal  from  this  detector  goes  to 
a  comparator  where  it  is  compared  with  the  average  intensity  which  has  been 
computed  from  the  electrical  signal  of  the  input  detector.  A  difference  in 
these  signals  causes  the  phase  shifter  driver  to  vary  the  nonreciprocal 
phase  shift  in  the  sensing  loop  until  the  output  and  average  intensities 
are  equal.  The  gyroscope  is  thus  maintained  at  the  stable  operating  point. 

,  The  signal  which  the  driver  applies  to  the  variable  non-reciprocal 
phase  shifter  is  measured  and  the  rotation  rate  is  computed.  Since  the 
gyroscope  is  at  or  near  a  stable  operating  point,  it  is  relatively  insen¬ 
sitive  to  reciprocal  birefringent  changes  in  the  sensing  loop. 
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The  polarization  controller  controls  the  slope 
response  curve  to  keep  it  from  going  to  zero  {environmental  fading).  The 
control  is  exerted  over  the  birefringent  parameters  9  and  0  in  Eq.  (8). 
The  degree  of  control  required  is  quite  loose  since  the  slope  does  not  need 
to  be  held  constant  at  a  particular  value,  but  just  needs  to  be  kept  away 
from  zero.  The  polarization  control  also  relaxes  the  limit  on  the  residual 
degree  of  polarization  of  the  source.  If  the  source  is  not  completely  un¬ 
polarized,  then  a  phase  measurement  error  appears  which  is  proportional  to 
the  degree  of  polarization.  This  phase  measurement  error  is  reduced  signi¬ 
ficantly  by  the  polarization  controller. 

All  components  except  for  the  variable  nonreciprocal  phase  shifter 
are  on  hand.  The  phase  shifter  still  needs  development,  before  it  can  be 
used  in  the  proposed  unpolarized  gyroscope. 

Advantages  of  this  unpolarized  gyroscope  are  reduced  noise  due  to 
coherent  backscattering  and  reciprocal  birefringence  changes,  linear  output 
over  a  broad  dynamic  range,  a  reduced  requirement  on  degree  of  polarization 
control,  and  omission  of  the  high  quality  polarizer  needed  for  polarized 
gyros. 

D.  NOISE  REDUCTION 

At  the  beginning  of  the  present  reporting  period  all  known  fiber  rota¬ 
tion  sensors  appeared  to  be  limited  in  their  minimum  measurable  rotation 
rate,  to  values  much  higher  than  the  theoretical  limit  due  to  quantum  noise. 
During  the  reporting  period  we  studied  this  problem  intensively  and  came  to 
the  conclusion  that  this  was  largely  due  to  coherent  Rayleigh  backscattering 
of  the  primary  circulating  waves.  This  is  scattering  from  the  material 
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inhomogeneities  in  the  fiber.  Such  inhomogeneities  are  characteristic 
of  all  fibers  now  fabricated,  and  this  phenomena  is  responsible  for  the 
principal  optical  propagation  loss  in  fibers.  However,  some  portion  of 
the  light  is  scattered  into  a  direction  so  that  it  is  still  confined  to 
the  core  and  propagates  along  the  fiber  in  the  opposite  direction  to  the 
original  signal.  This  occurs  in  all  fibers.  In  the  ring  Sagnac  sensor 
this  backward  scattered  signal  from  each  wave,  of  course,  propagates  along 
with  the  injected  opposite  traveling  wave,  is  guided  to  the  detector  where 
it  contributes  to  the  overall  noise  and  thus  degrades  the  signal-to-noise- 
ratio.  The  primary  limitation  to  rotation  sensitivity,  however,  is  the 
phase  error  produced  by  this  coherent  backscattering.  To  the  extent  that 
the  backscattered  waves  are  coherent  with  the  primary  signal  waves  they 
add  vectorially  at  unpredictable  phases  that  are  not  identical  for  the  two 
oppositely  directed  waves.  The  resulting  measured  phase  shift  is  indis¬ 
tinguishable  from  that  caused  by  rotation.  A  detailed  analysis  of  this 
scattering  and  its  effect  is  given  in  Appendix  D  which  is  the  text  of  a 
published  paper  describing  our  work  on  this  problem  done  under  this  contract. 

We  have  been  considering  various  methods  for  overcoming  this  problem 
and  describe  them  briefly  below.  Further  details  on  these  proposed  methods 
will  also  be  found  in  Appendix  D. 

The  phase  error  we  have  described  above  can  be  minimized  by  reducing 

the  coherence  between  primary  and  backscattered  waves.  One  method  of  ran- 

6“8 

domization  is  by  phase  modulation  of  the  loop.  If  the  phase  modulation 

^K.  Bbhm,  P.  Russer,  E.  Weidel,  and  R.  Ulrich,  Optics  Letters  64  (1981). 

^R.  A.  Bergh,  H.  C.  Lefevre,  and  H.  J.  Shaw,  Optics  Letters  6,  198  (1981). 

®R.  A.  Bergh,  H.  C.  Lefevre,  and  H.  J.  Shaw,  lOOC  '81  Technical  Digest,  p. 

116  (1981). 
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is  rapid  enough  to  allow  sampling  of  the  signal  at  a  frequency  M  ,  the 

-1/2 

cumulative  rotation  error  is  reduced  in  proportion  to  M  '  . 

A  greater  reduction  in  coherence  can  be  achieved  by  using  a  source 
having  a  short  coherence  length  (e.g.,  a  laser  diode).  ’  In  this  case 
only  those  scatterers  located  such  that  the  round  trip  transit  time  from 
the  input  differs  from  the  loop  transit  time  by  less  than  the  coherence 
time  of  the  source  will  make  large  contributions  to  the  phase  error.  This 
means  that  most  of  the  phase  error  is  due  to  those  scatterers  located  half¬ 
way  along  the  loop  within  a  fiber  segment  equal  in  length  to  the  coherence 
length  of  the  source.  Scatterers  located  in  other  parts  of  the  fiber  con¬ 
tribute  to  the  phase  error  to  a  lesser  degree,  the  extent  of  which  is  deter¬ 
mined  by  the  coherence  properties  of  the  source. 

These  methods  have  been  tried  in  another  program  at  Stanford  and  have 
increased  fiber  gyro  sensitivity  by  orders  of  magnitude,  but  significant 
error  still  remains.  One  way  to  further  reduce  this  error  is  to  use  signal 
pulses  that  are  short  compared  to  the  loop  transit  time.  (Alternatively, 
one  could  use  longer  pulses  that  are  chirped  to  form  the  spectral  equivalent 
of  the  short  pulses,  thereby  avoiding  peak  power  problems.)  The  output  of 
the  loop  equal  to  the  pulse  length  in  the  fiber  contributes  to  the  error. 

While  the  use  of  short  pulses  may  appear  to  yield  the  same  error  reduc¬ 
tion  as  would  the  use  of  a  source  having  an  equally  short  coherence  time, 
pulsed  operation  in  fact  has  several  advantages.  The  offending  segment  of 
the  loop  is  more  clearly  defined  by  a  signal  pulse  than  by  the  source  co¬ 
herence  time.  Furthermore,  the  self-gating  aspect  of  pulse  operation 

q 

K.  Bbhm,  P.  Marten,  K.  Peterman,  and  E.  Weidel,  lOOC  ‘81,  Post  Deadline 
Paper  TuL7  (1981 ) . 

A.  Bergh,  H.  C.  Lefevre,  and  H.  J.  Shaw  (submitted  for  publication). 
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insures  that  scatterers  located  elsewhere  in  the  loop  make  no  contribution 
to  either  the  phase  error  or  the  shot  noise.  This  allows  us  to  concentrate 
on  reducing  the  backscattering  efficiency  of  only  this  short  segment  of 
fiber  as  a  means  of  reducing  the  backscattering  error  even  further  (Fig. 
10).  We  propose  that  this  segment  of  the  loop  be: 

(1)  index  matched  to  reduce  the  capture  efficiency  of  the  fiber, 

(2)  replaced  by  an  ultra-pure  fiber  segment  having  few  scatterers, 

(3)  replaced  by  a  hollow  light  guide,  or 

(4)  evacuated  entirely  to  insure  that  no  sources  of  scattering  or 
reflection  are  present. 

Finally,  our  analysis  indicates  that,  all  other  parameters  being 
constant,  the  use  of  a  longer  signal  wavelength  reduces  the  phase  error 
due  to  backscattering  by  roughly  (Fig.  11). 

It  is  probably  that  by  implementing  these  techniques  the  sensitivity 
of  existing  fiber  gyros  can  be  substantially  improved. 
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MODIFIED  FIBER  SEGMENT 


FIGURE  10 


a  (X)  /a  ( Xo) 


FIGURE  n 
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I.  INTRODUCTION 

Efforts  are  being  made  in  various  laboratories^  to  develop  fiber 
optic  rotation  sensors  (FORS)  as  alternatives  to  mechanical  gyroscopes  and 
ring  laser  gyros  (RLG)  in  some  applications.  All  FORS’s  utilize  the  Sagnac 
effect^^  in  a  ring  interferometer.  The  ring  interferometer  is  formed  from 
a  length  of  "single"  mode  optical  fiber  wound  N  times  around  an  area  A  . 

In  the  presence  of  rotation,  countra-directional  optical  waves  propagating 
around  the  ring  interferometer  experience  a  phase  shift  (Sagnac  Phase  shift 
(J)^)  given  as 


Stt  NA 
Xc 


(1) 


where  X  and  c  are  the  free  space  values  of  wavelength  and  speed  of  light 
and  is  the  rotational  velocity  of  the  ring  interferometer.  In  most  FOR's 

the  intensity  of  the  interferometer  pattern  of  the  ring  interferometer  is 
monitored  in  a  localized  area.  Changes  in  the  rotation  rate  of  the  FORS 
give  rise  to  change  in  the  Sagnac  phase  shift  which  in  turn  is  detected  as 
a  change  in  intensity  in  the  monitored  area  of  the  interference  pattern. 

While  the  fringe  pattern  appears  relatively  stable  to  the  eye,  it 
actually  possesses  variations  with  time  which  are  sources  of  error.  Environ¬ 
mental  effects  otlicr  than  rotation,  such  as  temperature,  vibration,  etc., 
can  cause  variations  in  both  the  positions  and  intensities  of  the  fringe 
pattern,  wliich  often  cannot  be  distinguished  from  the  desired  changes  due 
to  rotation.  It  has  been  shown  ’  that  by  using  polarized  light  at  the 
input  to  the  interferometer,  and  using  polarization  filters  in  the  output, 
the  operation  of  the  sensing  loop  can  be  made  completely  reciprocal,  and 
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the  positions  of  Che  interference  fringes  are  then  stable  against  environ¬ 
mental  changes  which  are  homogeneous  over  the  space  occupied  by  the  sensing 

loop.  However,  the  intensity  and  contrast  of  the  fringes  can  still  be 

13  14 

environmentally  affected.  Work,  is  being  done  ’  on  approaches  to  reducing 
these  effects. 

In  the  following  sections  we  examine  further  the  effects  of  polariza¬ 
tion  on  the  operation  of  a  FORS.  This  includes  theoretical  and  experimental 
investigation  of  operation  using  polarized,  unpolarized,  and  partially 
polarized  light. 

II .  THEORY 

Real  "single"  mode  fibers  propagate  two  orthogonal  polarization  modes. 
These  two  modes  have  sliglitly  different  propagation  velocities.  Environ¬ 
mental  effects  mix  the  two  polarization  modes.  This  mixing  of  the  ortho¬ 
gonal  modes  can  result  in  the  instabilities  observed  in  FORS's  if  proper 
care  is  not  taken.  It  is  important  for  the  theoretical  analysis  to  realize 
that  while  the  state  of  polarization  can  change  in  a  real  '"single"  mode  fiber 
the  degree  of  poJ.arization  is  unchanged.  This  allows  us  to  treat  the  fiber 
as  a  four  port  linear  system.  In  Fig.  1,  the  four  ports  are. shown  explicitly 
Although  the  polarization  states  are  labelled  x  and  y  they  can,  in 
general,  be  any  two  orthogonal  polarization  states.  Also  shown  schematically 
in  Fig.  1  are  the  possible  channels  between  dlie  four  ports.  A  cliannei  is 
characrcrized  by  a  complex  scattering  coefficient.  The  matrices  describing 
the  scattering  between  the  fiber  ends  among  the  polarization  modes  are 
recognized  as  Jones  matrices  for  a  birefringent  element.  This  representation 
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End  1 


End  2 


Ports 


Fig.  1-Four-port  linear  system  model  of  a  single  mode  fiber. 
The  ports  are  x-j ,  -^1  *  -^2  •  channels  between 

the  ports  are  shown  schematical ly  along  with  their 
complex  scattering  coefficients. 
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of  a  real  "single"  mode  fiber  by  Jones  matrices  was  first  presented  by 
Ulrich,  et  al.^^ 

In  Fig.  2,  the  system  which  we  will  analyze  is  shown.  We  have  a 
length  of  single  mode  fiber  wound  N  times  around  an  area  A  .  We  repre¬ 
sent  the  beam  splitting  and  recombining  elements  as  a  black  box.  We  will 
neglect  the  details  of  the  beam  splitting  and  recombining  in  our  calculations. 
We  should  point  out,  however,  that  these  details  can,  in  principle,  be  a 
source  of  nonreciprocity  in  the  system  unless  properly  configured.  Inclusion 
of  these  details  can  give  rise  to  scale  factor  changes,  a  phase  shift,  reduc¬ 
tion  in  contrast,  and  a  change  in  the  degree  of  polarization.  We  are 
interested  here  only  in  the  effects  of  the  fiberoptic  sensing  loop  on  the 
system.  We  assume  that  the  laser  can  be  operated  with  any  degree  of  polari¬ 
zation  required.  In  Fig.  2  a  general  polarization  analyzer  can  be  inserted 
before  the  detector  to  pass  a  given  polarization  state  only. 

The  Jones  matrix  for  propagation  from  end  1  to  end  2  (Fig.  1)  is  , 

and  for  propagation  from  end  2  to  end  1  is  •  This  Jones  formalism 

allows  the  inclusion  of  the  birefringent  effects  which  are  in  general  present 
in  real  single  mode  fibers.  The  matrices  and  T^2  unitary  for 

zero  rotation  rate.  With  rotation  present  the  explicit  forms  for  T2j^  and 


"  t(x2,Xj^) 

e 

t(y2.5<i) 

t(y2.y2^) 

-j't'  /2 

^tCx^.x^) 

t(.\^,y2)' 

s 

e 

t(yj^>y2)- 

(2) 


(3) 


4 


2--Schematic  of  system  analyzed. 


Here  t})  is  the  Sagnac  phase  shift  as  given  by  Eq.  (1).  The  t's 
s 

(Fig.  1)  are  the  complex  scattering  coefficients  of  the  fiber,  and  are 
environmentally  dependent.  The  principle  of  reciprocity  imposes  the 

4' 

following  condition 

t(k,ii)  =  t(Jl,k)  (4) 

This  means  that  the  effect  of  a  channel  is  independent  of  the  direction  in 
which  the  channel  is  traversed.  Equation  (4)  is  valid  only  if  there  is  no 
rapidly  changing  gradient  along  the  fiber.  By  rapid  we  mean  fast  compared 
to  the  transit  time  of  the  fiber.  In  our  calculations,  we  will  assume  that 
the  fiber  is  lossless  and  that  the  optical  power  propagating  in  each  direc¬ 
tion  in  the  fiber  are  equal. 

The  input  into  each  end  of  the  fiber  is  represented  as  a  vector  E^^  , 
where  the  components  of  the  vector  are  the  E  fields  of  the  two  polariza¬ 
tion  modes.  Allowing  the  input  to  propagate  around  the  fiber  in  opposite 
directions  and  recombining  the  outputs  from  each  end  of  the  fiber,  we  obtain 


E 

out 


E. 

in 


+  T 
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E. 

in 


(5) 


This  signal  can  be  passed  through  a  general  polarization  analyzer 
represented  by  a  matrix  P  .  The  detected  E  field  is  then 


-*■ 

E,  =  P  E 
det  out 


(6) 


We  will  now  calculate  +  det 
analyzers . 


E,  for  various  E.  and  polarization 

det'  in 
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III.  CALCULATED  CASES 


Case  1.  We  will  let  E.  be  polarized  in  the  first  state  of  polari- 


in 


zation.  The  polarization  analyzer  will  pass  this  state.  The  forms  for 


E .  and  P  are 
in 


E.  =  {1,0} 

in 


(7) 


P  = 


1  0 
0  0 


(8) 


We  now  perform  the  operation  in  Eqs.  (5)  and  (6)  and  calculate  the 


normalized  detected  intensity  I 


det 


=  l/2l  tCx^.x^)  1^(1  +  cos  c|)^)  (9) 

Figure  3  shows  a  plot  of  I,  versus  4  for  various  environmental 
conditions,  i.e.,  values  of  t(x^,x^)  .  The  environment  conditions  are 
parameterized  by  9  and  ({)  .  9  is  a  measure  of  the  power  transfer  between 

the  two  polarization  modes.  <f)  is  also  the  non-reciptocal  phase  shift 
betv.’een  the  contradirectional  waves  which  enter  the  fiber  in  one  polarization 
and  exit  the  fiber  in  the  other  polarization.  Refer  to  Appendix  A  for  the 
derivation  6  and  <j)  as  environmental  parameters.  The  plots  are  sinusoidal. 
The  minimum  values  of  Che  plots  (i.e.,  I,  =0)  are  insensitive  to  enviroii- 
mental  changes.  Unfortunately  the  sensitivity  to  rotation  at  the  minima  is 
also  zero.  Tlie  maxima  of  the  plots  are  environmentally  sensitive.  The 
value  of  the  maximum  varies  between  0  and  1.  The  position  (i.e.,  the  value 
of  which  maximum)  is  not  environmentally  dependent.  This 

means  tliat  various  detection  schemes  (i.e.,  detection  in  quadrature)  can  be 
used  to  compensate  for  die.  environmentally  produced  intensity  changes. 
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Normal i zed 
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-Normalized  intensity  vs  Sagnac  phase  shift  for  polarized  input 
Analyzer  set  to  pass  input  polarization.  Parameter  6  and  ij) 
are  discussed  in  the  text  below  Eq .  (9). 


These  schemes  give  large  errors  when  tCx^.x^)  0  .  To  prevent  this 
polarization  conserving  fiber  or  a  means  of  active  polarization  control  in 
the  fiber  is  necessary. 

Case  2.  We  will  let  E.  be  defined  as  in  Case  1,  but  the  polariza- 
-  in 

tion  analyzer  will  now  pass  the  orthogonal  polarization  state.  The  form  for 
P  is  now 


P 


0  0 
0  1 


(10) 


Using  Eqs.  (7)  and  (10)  in  Eqs.  (5)  and  (6)  we  obtain  the  normalized 
detected  intensity 


=  1/4|  t(y2,Xj^)  !  ^  +  l/4l  c(yj^,X2)  1^  +  l/2l  t(y2,x^)  I  1  t(y^,X2)  1 

(11) 

X  cos{<!)^  +  arg[t(y2,x^)  ]  -  arg[t(Xj^,X2)  ] } 

Figure  4  shows  plots  of  versus  for  the  same  environmental 

conditions  as  was  used  for  Fig.  3.  All  plots  in  this  paper  will  use  the 
same  set  of  environmental  conditions.  See  Appendix  A  for  the  analytical 
form  of  the  T's  used  in  the  calculations.  The  plots  are  again  sinusoidal. 
The  position  of  the  minima  and  maxima  are  now  environmentally  sensi¬ 
tive.  Tlic  value  of  the  minima  is  still  environmentally  insensitive. 

The  value  of  the  maximum  is  environmentally  sensitive.  Because  the  posi¬ 
tions  of  the  minima  and  maxima  can  change,  there  is  no  detection  scheme 

which  will  allow  one  to  measure  tf)  .  This  case  is  an  incorrect  mode  of 

s 

operation  of  a  FORS, 
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,  except  that  analyzer  is  set  to  reject  input  polarization. 


Case  3.  E 


is  the  same  as  in  the  two  previous  cases.  There  is  no 


in 

polarization  analyzer  before  the  detector.  In  this  case  the  normalized 
detected  intensity  is  computed  to  be 


Ijgj.  =  l/2lt(x^,x^)I^(l  +  cos  +  l/4|t(x^,x^)|^  +  l/4lt(y^,X2)|^ 

+  1/2|  tCy^.x^)  I  I  tCy^.x^)  i  cos{(|)^  +  arg[t(y2,x^)  ]  -  arg[t (y^,X2)  ] } 

(12) 

Figure  5  shows  plots  of  I,  versus  6  for  the  same  set  of  environ- 

det  s 

mental  conditions  as  before.  The  positions  of  the  ;miniiiia  and  maxima 
vary  with  the  environment.  This  case  is  also  an  incorrect  mode  of  operation 
of  a  FORS .  It  is  interesting  to  note  that  in  Fig.  5  the  average  value  or 
"DC"  level  of  the  plots  are  constant  at  a  value  of  1/2.  The  minima  of 
the  plots  are  now  environmentally  sensitive. 

Case  4 .  We  will  now  treat  the  case  where  the  input  light  is  com¬ 
pletely  unpolarized.  We  will  represent  the  input  vector  as 


E. 

in 


E  (t)e 

X 


j'^‘(t)/2 


E  (t)e 

y 


-j4)(t)/2 


(13) 


(P  is  a  random  variable  uniformly  distributed  over  the  interval  (-it, it)  . 

E  (t)  and  E  (t)  are  the  instantaneous  fields  of  the  two  polarization 
X  y 

components.  The  time  averaged  values  of  tJie  squared  magnitudes  of 
and  E  (t)  are 

y 

|E^Ct)I^  =  |EyCt)I^  =  1/2  (14) 


v.'liere  I  is  tlie  total  intensity  of  tlic  light.  There  is  no  polarization 
analyzer  before  the  detector.  Substituting  Eq .  (12)  into  Eq.  (5),  finding 
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the  magnitude  squared  of  E 


,  taking  the  time  average,  and  performing 


det 

some  algebraic  manipulations  we  obtain  for  the  normalized,  time  averaged, 
detected  intensity 

Ijgj.  =  1/2  +  1/2(1  -  I/2I  t(x2,yj^)  -  tCx^.y^)  l^)cos  (15) 

This  equation  is  valid  if  the  coherence  time  of  the  polarization  of  the 

light  source  is  much  shorter  than  the  response  time  of  the  detection  system 

used.  Figure  6  shows  plots  of  I,  versus  c()  for  the  standard  set  of 

det  s 

environmental  conditions  used  before.  The  position  of  the  minima  and 
maxima  are  again  environmentally  insensitive.  The  average  value  of  the 
detected  intensity  is  again  constant  at  a  value-  of  1/2  . 

The  most  important  point  to  realize  is  the  behavior  of  the  detected 
intensity  at  the  77/2  points  (i.e.,  “  ro  ^/2,  ra  =  +  1,  +  3,  ...).  If 

the  FORS  is  biased  to  operate  at  a  Tr/2  point,  there  will  be  _r^  environ¬ 
mental  dependence  of  the  detected  Intensity.  The  FORS  will  be  stable.  The 
sensitivity  or  scale  factor  at  the  Tr/2  point  is  still  environmentally 
dependent.  This  environmental  dependence  can  be  handled  in  various  ways 
which  are  outside  the  scope  of  this  paper. 

In  practice  it  is  difficult  to  obtain  sources  of  completely  unpolarized 
light  of  sufficient  brightness  and  of  such  a  geometry  as  to  couple  efficiently 
into  a  single  mode  optical  fiber.  This  problem  can  be  overcome  by  using  two 
linearly  polarized  sources  which  are  uncorrclated.  The  plane  of  polarization 
of  Llie  sources  arc  at  right  angles  to  each  other.  Anotiier  method  is  to  use 
a  linearly  polarized  source  wliose  plane  of  polarization  rotates  rapidly.  A 
mucli  more  serious  problem  is  that  tlie  optical  elements  exterior  to  the  loop 
arc  in  genera]  polar ization  sensitive.  These  olemeiit.s  tend  to  partially 
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Normal ized 
Intensity 


lA 


-Same  as  Fig.  3,  except  that  input  is  completely  unpolarized,  and  no 
analyzer  is  used. 


polarize  the  input  light.  Because  of  these  two  problems  in  obtaining 
completely  unpolarized  light  at  the  input  to  the  fiber  optic  sensing  loop, 
it  is  necessary  to  determine  the  behavior  of  the  FORS  for  partially  polarized 
input . 

Case  5.  Partially  polarized  input  with  degree  of  polarization  P  . 

For  simplicity  we  assume  that  the  degree  of  polarization  P  of  the  contra- 
propagating  waves  are  equal. A  wave  of  degree  of  polarization  P  can 
be  decomposed  on  an  intensity  basis  into  a  completely  polarized  wave  (P  =  1) 
and  a  completely  unpolarized  wave  (P  =  0) .  If  total  intensity 

of  the  partially  polarized  wave,  then  the  intensities  of  the  polarized  and 
unpolarized  components  of  the  wave  are 


I  1  -  .  =  P  I  ^ 

polarized  cot 


I  1  •  j  . 

unpolarized  tot 


(16) 

(17) 


The  response  of  the  FORS  to  a  partially  polarized  input  can  be  expressed 
as  the  sum  of  the  responses  of  the  FORS  to  the  polarized  and  unpolarized 
components  of  the  partially  polarized  input. 

Let  there  be  no  analyzer  before  the  detector.  The  response  is  obtained 
by  combining  Eqs.  (11),  (15),  (16)  and  (17).  The  normalized  detected 
intensity  can  be  expressed  as 


I 


det 


I  ZltCx^.x^)!^  cos^  tCy^.x^)  I  ^  +  |j  t(yj^,X2)  j  ^ 

P  1 

+  -jl  t(y2,Xj^)  I  1  t(yi,X2)  Icos  +  ar2[t(y2,x^)  ]  -  ar  g  [  t  ( ,  x^)  ] 

+  1  +  (1  -  I  t(x  ,x  )  -  t(x  ,x J  (^)cos  4)  I  (18) 

2  '21  12'  s| 
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Figures  7,  8,  and  9  are  the  plots  of  I,  versus  i  from  Eq,  (18) 

det  s 

for  P  =  0.25  ,  0.50  ,  and  0.75  respectively.  The  same  environmental  con¬ 
ditions  are  used  as  before.  The  average  value  of  the  detected  intensity 
is  still  constant.  The  position  of  the  maxima  and  minima  are  environ¬ 
mentally  sensitive.  The  degree  of  sensitivity  increases  as  the  degree  of 
polarization  increases.  For  partially  polarized  inputs  the  7t/2  points 
are  no  longer  stable.  We  now  ask  what  maximum  degree  of  polarization  is 
allowable  for  a  given  sensitivity  to  rotation  of  a  FORS,  assuming  that  the 
environmentally  produced  phase  and  amplitude  shifts  in  the  sensing  loop  are 
the  limiting  factors.  We  consider  a  FORS  biased  to  the  Tr/2  point.  Assuming 
no  noise,  and  a  degree  of  polarization  P  =  0  ,  we  rotate  the  gyroscope  at 
some  angular  velocity  Q  ,  and  calculate  the  intensity  change  I(n)  .  We 
next  let  the  degree  of  polarization  be  nonzero.  We  then  calculate  the  maximum 

environmentally  produced  intensity  change  I  (Environment)  at  zero  rotation 

max 

rate  (f2  =  0)  .  We  then  ask  what  is  Che  value  of  the  degree  of  polarization 

P  for  v;hich  I(oQ)  =  I  (Environment).  This  value  of  P  is  the  maximum 

max 

value  of  P  allowable  for  a  minimum  rotation  sensitivity  ^  .  We  assume  a 

2 

fiber  optic  sensing  loop  enclosing  an  area  of  1  m  with  250  turns  and  operated 

O 

at  6323  A  .  We  assume  a  desired  sensitivity  of  3  deg/hr,  which  is  about  an 

order  of  magnitude  better  than  has  been  obtained.  This  rotation  rate  corres- 

“3 

ponds  to  a  Sagnac  phase  shift,  Eq .  (1),  of  0.3  x  10  rad  .  For  these 
conditions  we  obtain  a  maximum  degree  of  polarization  of  8  x  10  ^  .  This  is 
quite  small,  and  may  be  difficult  to  obtain  in  practice.  If  one  is  able  to 
control  the  state  of  the  polarized  part  of  Che  partially  polarized  wave  in 
the  fiber  optic  sensing  loop,  the  limit  on  degree  of  polarization  is  reduced. 
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-Same  as  Fig.  3  except  that  input  is  partially  polarized  with  degree  of 
polarization  P  =  0.25  and  no  analyzer  is  used. 
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-Same  as  Fig.  7  except  p  =  0.50 


Norma 

Inte 
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-Same  as  Fig.  7  exc 


Assume  that  the  fiber  maintains  linear  polarization  but  rotates  the  plane 
of  polarization.  Under  this  condition  if  the  plane  of  polarization  can  be 
stabilized  to  ^  2.5  degrees,  the  maximum  degree  of  polarization  becomes 
0,05  (5%)  for  a  3  deg/hr  sensitivity. 

The  appendix  contains  explicit  forms  for  the  Jones  matrices.  It  also 
contains  a  new  formalism  for  the  polarization  analysis  of  FORS's  using  the 
coherency  matrix.  This  formalism  places  the  polarization  analysis  of  FORS's 
on  a  firm  theoretical  foundation,  and  is  especially  useful  for  the  case  of 
partial  polarization. 

IV.  EXPERIMENTS 

Experiments  have  been  performed  using  a  CW  fiber  gyro  with  a  single 
mode  fiber  of  600  m  length  and  1  ra  diameter,  excited  with  a  helium  neon 
laser.  Detector  output  current  was  displayed  on  an  oscilloscope  as  a  function 
of  rotation  rate  of  the  system,  giving  patterns  of  the  types  plotted  in  Figs. 

3  to  9.  Environmental  effects  were  simulated  by  subjecting  a  short  length 
of  fiber,  inside  the  loop,  to  a  double  twist.  Specifically,  Che  fiber  was 
clamped  Co  the  rotating  table  at  two  points  spaced  about  18"  along  the  fiber. 
The  fiber  segment  between  the  clamps  contained  slack,  and  a  rotatable  clamp 
was  attached  to  Che  center  of  this  segment,  which  could  be  rotated  about  the 
fiber  axis.  The  fiber  was  twisted  through  an  angle  which  produced  approxi¬ 
mately  LliG  maximum  change  in  the  displaced  output  pattern.  This  angle  was 
generally  in  the  range  of  90°  to  180°.  The  actual  value  was  not  important 
as  the  purpose  was  to  confirm  the  qualitative  predictions  of  the  theory, 
wliicli  sliow  some  polarization  states  to  be  free  of  certain  homogeneous 
environmental  effects,  wliile  others  are  not. 
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Figure  10  is  for  the  case  of  linear  polarized  input  light  together 
with  a  linear  polarization  filter  at  the  detector,  set  to  pass  only  the 
same  polarization.  The  lower  trace  corresponds  to  zero  twist  angle  and 
the  upper  trace  results  from  introducing  the  twist-  Note  that  the  positions 
of  the  maxima  and  minima  do  not  shift  along  the  horizontal  axis,  confirming 
that  recpirocity  is  satisfied.  Note  also  that  the  peak  amplitude  does 
change,  but  that  the  amplitude  of  the  minima  does  not,  as  predicted. 

Figure  11  shows  the  results  when  the  polarizer  is  rotated  to  pass 
only  the  polarization  orthogonal  to  the  input  polarization,  and  Fig.  12 
results  when  the  polarizer  is  removed.  Lovzer  and  upper  traces  are  for  the 
same  tv/ist  states  as  before.  Note  in  both  cases  the  positions  of  maxima 
and  minima  shift  horizontally  as  a  function  of  twist,  i.e.,  the  system 
is  nonreciprocal,  and  amplitude  variation  with  twist  is  again  present,  all 
as  predicted. 

Figure  13  results  when  unpolarizcd  light  is  used,  and  no  polarizers 
are  employed.  The  upper  trace  corresponds  to  zero  twist,  while  the  lower 
traces  result  when  twist  is  introduced.  Note  that  (1)  the  horizontal  posi¬ 
tions  of  the  maxima  and  minima  remain  fixed  and  (2)  the  positions  of  the 
tt/2  points  remain  fixed,  again  as  predicted  by  the  theory. 
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Fig.  10 — Oscillogram  of  detector  current  (vertical  axis)  vs 
rotation  rate  (horizontal  axis) ,  for  linear  input 
polarization  and  co-polarized  detector  filter. 
Lower  trace:  fiber  twist  absent.  Upper  trace: 
fiber  twist  present. 


I'lg.  11 — Same  as  Fig.  10,  except  detector  filter  was  cross- 
polarized  . 


Fig-  12 — Same  as  Fig.  10  except  no  detector  polarization  filtel" 
was  used. 


Fig.  13 — Oscillogram  of  detectoic  cui-ronu  (vertical  axis)  vs 
rotation  rate  (lior  Xzontal  nxi;'),  for  unpo  J  arizeci 
input  and  no  polarization  rj:._cr;'.  op.  cr  trace; 
fiber  twist  ab.scnt.  Lower  traces.;  fiber  twist 
prosen  f , 
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V.  APPENDIX 


We  now  consider  first  an  explicit  form  for  the  Jones  matrix  for  a 
real  single  mode  fiber.  Pictured  schematically  in  Fig.  A-1  is  a  real  single 
mode  fiber  which  propagates  two  polarization  modes  of  propagation  constants 


*1 

and  respectively.  At 

a  point 

z  along  the 

fiber  there 

exists  a  . 

linear  scattering  center.  This 

scattering  center  can 

exchange  all  or 

part 

of 

the  powers  in  the  two  modes. 

At  this 

point 

the  Jones  matrix 

T 

BA 

can 

be 

represented  as 

cos  0 

^'*’11 

sin  9 

T  = 

BA 

e 

e 

(19^ 

-sin  6 

cos  0 

3*22 

e 

e  > 

The  Jones  matrix  T,„ 

AB 

can 

be 

obtained  from  T^^  by  using 

the  principle  of 

reciprocity . 

It  now  remains 

to  specify  the  forms  of  the  6 . . ' s  . 

3-J 

Here  (t .  .  reprs' 
ij 

sents  the  phase  shift 

imparted  to  a  wave  which  started  in 

the  polariza- 

tion  and  finished  in 

the  i 

th 

polarization.  By  inspection 

of  Fig.  1  we  can 

write 

(a) 

hi 

= 

61L 

(b) 

'^22 

= 

(c) 

^12 

-A  3  2  +  3j|^L 

(20) 

(cl) 

^21 

= 

A3^  +  32^ 

(c) 

A0 

= 

3.  -  32 

(f; 

S 

-■ 

(Bj_  +  0^)/2 
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By  redefining  the  origin  of  phase  (i.e.,  $  if  “  3  L)  we  can  write 


(19)  as 


r  i  ^  L 

J  2  ^ 

cos  9  e 

sin  6  e 

T 

.A3/L  \ 

.  A6  , 

BA 

.  e  -J  (i-V 

sin  e 

2  ^ 

cos  6  e 

(21) 


Finally  we  call  A31.  =  and  A3  =  *fpQ/2  (PD  =  path  difference). 

The  explicit  forms  for  the  therefore 


(a) 

^(x^, 

.x^) 

=  cos 

0 

(b) 

t(x^, 

'V 

=  sin 

9 

e 

j<i> 

(c) 

tCy^. 

Xi) 

=  -sin 

0 

e 

(d) 

t(y2. 

.yp 

=  cos 

0 

(22) 


For  a  fiber  with  N  scattering  centers  the  Jones  matrix  still  has 
elements  of  the  above  form. 

We  next  develop  an  alternative  formulation  of  the  effect  of  the 
optical  fiber  loop  on  FORS  performance.  Consider  a  wave  of  arbitrary  polari¬ 
zation  propagating  in  the  z-direction,  given  by 


E(z,t)  =  I 


(23) 


The  coherency  matrix  (J)  is  defined  as 


<rj->  = 


<\<>  <Vy> 


L 


i:.  > 

\  V  r.  / 


(2/0 
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time  average 


+  ;  hermit ian  conjugate 
E  E  :  dyadic  product 
*  :  complex  conjugate 

The  intensity  (I)  of  the  wave  represented  by  J  is  given  as 


I  =  tr  J 

The  off-diagonal  elements  of  J  are  the  cross-correlation  of  the  two 
polarization  components.  J  is  hermitian. 

Examples  of  the  coherency  matrix  for  various  polarized  waves  of 
unit  intensity  are  as  follows: 


(a) 

Linear  polarization 

'l  0 

in  the  x  direction 

J 

_0  0_ 

(b) 

Linear  polarization 

'l/2 

111 

at  45*’  to  x-axis 

J 

_l/2 

1/2. 

(c) 

Right-hand  circular 

■ 

1/2 

j/2‘ 

polarization 

J 

-j/2 

1/2  _ 

'  1/2 

0 

(d) 

unpolarized  light 

J 

_  0 

1/2  , 

(c) 

50%  partial  polarized 

light-polarized  part 

J 

3/4 

0 

along  x-nxis 

-  0 

/14  _ 

(25) 


For  propagation  of  the  coherency  matrix  through  a  linear  systems,  let  the 
system  be  described  by  the  matrix  L  ,  where  L  is  defined  by  the  relation 


E  =  L  E. 
out  in 


(27) 


The  coherency  matrix  J  is  found  from  the  coherency  matrix 

out 

J.  by  the  following  expression 


J 


out 


(28) 


For  a  ring  interferometer  operated  as  a  Sagnac  rotation  sensor,  we  can 
express  L  as  (see  Fig.  11): 


reciprocal  BS  21  BS  BS  12  BS 


(29) 


L  .  T  =  P(Rr^T_R^_  +  T„„T„T  ) 

nonreciprocal  BS  21  BS  BS  12  BS 


(30) 


P  ,  defined  before  and 


BS 


fr  O' 

1 

rt 

O 

X 

T  = 

X 

0  r 

^BS 

0  t 

L  y  J 

L  y  J 

(31) 


v;herc  r  ,  r  ,  t  ,  t  are  complex  valued  transfer  coefficients  for  the 
X  y  X  y 

E-field. 


To  duplicate  the  results  obtained  earlier,  set 
(the  identity  matrix).  For  cases  1,  2,  and  3  use  ”  0  o] 


use  J .  = 

in 


'l  o' 

and  for  case  5  use 

Y  (1-P) 

1 

0 

.0  1. 

0 

1 

2 

(1  +  P)  - 

For  case  4 


where  P 


is  tlie  degree  of  polarization  of  the  wave.. 
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APPENDIX  B 


COHERENCE  MATRIX  FORMULISM 

The  description  of  the  polarization  of  a  light  wave  requires  the 
specification  of  its  state  and  degree  of  polarization.  The  coherence 
matrix^  of  the  light  wave  contains  this  information.  Any  two  ortho¬ 
normal  polarization  states  can  be  used  as  the  basis  for  the  coherence 
matrix.  For  the  analysis  performed  in  this  appendix,  the  two  states 
are  chosen  as  linearly  polarized  waves  along  the  x  and  y  axes.  Let 
E^  and  be  the  instantaneous  projections  of  the  electric  field,  E  , 
of  the  light  wave  along  the  x  axis  and  y  axes  respectively.  The  co¬ 
herence  matrix  J  of  the  light  wave  is  given  in  Eq.  (1). 

/<Ex  <>  <^x 

J  =  (1) 

\<S  <>  Ol 

The  brackets  ^ >  signify  the  infinite  time  average  and  the  asterisks 
signify  the  complex  conjugate.  There  are  several  properties  of  the  coher¬ 
ence  matrix  which  are  useful.  The  first  property  is  that  the  intensity 
(I)  of  the  light  wave  is  the  trace  of  the  coherence  matrix.  Second,  the 
coherence  matrix  is  hermitian  and  third,  the  degree  of  polarization  (P) 
of  the  light  wave  can  be  calculated  from  the  eigenvalues  (X-j  ,  X2)  of 
the  coherence  matrix.  These  properties  are  summarized  in  Eq.  (2). 

I  =  TrJ 

Jki  = 

P  =  ix^  -  X2l/(X^  +  X^) 
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The  coherence  matrices  for  various  states  and  degrees  of  polarization  are 
shown  in  Fig.  1.  Linear  polarization  along  the  x  axis  and  the  y  axis 
form  the  basis  set  for  these  coherence  matrices.  The  off-diagonal  terms  of 
the  coherence  matrix  represent  the  correlation  between  the  fields  along  the 
basis  vectors. 

Suppose  that  we  have  a  black  box  as  shown  in  Fig.  2.  A  light  wave 
characterized  by  is  input  into  the  box.  The  output  light 

wave  from  the  box  is  characterized  by  =  ^^x’^y^out  ‘ 

tains  only  linear  components,  the  output  can  be  related  to  the  input  by  a 
2x2  matrix,  S  ,  of  complex  numbers.  The  mathematical  relation  is  given 
in  Eq.  (3). 


(3) 


If  the  input  and  output  fields  are  related  by  Eq.  (3),  it  is  easy  to  show 
that  the  input  and  output  coherence  matrices  are  related  by  the  formula 
in  Eq.  (4)  where  represents  the  hermitian  adjoint  of  S 


'^out 


(4) 


Equation  (4)  is  valid  only  if  the  matrix  S  describing  the  linear  system 
does  not  change  appreciably  during  the  interval  of  time  over  which  the 
time  averages  contained  in  J  are  taken.  It  is  also  easy  to  show  that  if 
two  light  waves  with  coherence  matrices  J-j  and  J2  are  superimposed,  the 
coherence  matrix  of  the  total  light  wave  is  just  the  sum  of  the  coherence 
matrices  of  each  light  wave. 

This  formalism  is  an  important  advance  because  it  can  be  used  to  treat 
quasi-monochromatic  light  (i.e.,  light  with  degree  of  polarization  other 
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Figure  1.  Coherence  matrices 
of  polarization. 


Linearly  polarized  along 
X  axis 


Linearly  polarized  along 
Y  axis 


Linearly  polarized  at  45 
to  X  axis 


Right  hand  circular 
polarization 


Left  hand  circular 
polarization 


Unpolarized  light 


Partially  polarized  light  with 
polarized  part  being  linearly 
polarized  along  X  axis.  Degree 
of  polarization  is  1/2. 


for  various  states  and  degrees 
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Figure  2.  A  linear  system  is  represented  as  a  black  box  by 
its  transfer  matrix,  S.  The  transfer  matrix  relates  the 
output  vector  to  the  input  vector. 
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than  unity).  Previous  formalisms  could  treat  only  completely  polarized 
light.  The  analysis  of  the  unpolarized  fiber  gyro  is  greatly  facilitated 
by  this  formalism. 
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APPENDIX  C 


JONES  MATRICES  FOR  BIREFRINGENT,  SINGLE  MODE  FIBER  GYROS 


The  birefringence  of  single  mode  fiber  gyros  is  described  by  Jones 
matrices.  This  appendix  develops  Jones  matrices  for  the  optical  fiber, 
the  polarizer,  and  the  beamsplitter.  The  Jones  matrix  of  the  fiber  gyro 
can  be  expressed  in  terms  of  the  Jones  matrices  of  the  components. 

Describing  an  optical  fiber  by  Jones  matrices  is  equivalent  to  model¬ 
ing  the  optical  fiber  as  a  general  four  port,  linear  system  as  shown  in 
Fig.  1  with  the  four  ports  and  the  paths  between  the  ports  depicted  expli¬ 
citly.  Each  path  is  characterized  by  a  complex  scattering  coefficient. 

The  complex  scattering  coefficient  describes  the  phase  shift  collected  by 
the  light  wave  and  the  amplitude  of  the  light  wave  which  travels  a  specific 
channel  through  the  system.  The  scattering  coefficients  are  dimensionless. 
We  can  write  a  Jones  matrix  S12  for  propagation  from  end  number  2  to  end 
number  1  of  the  fiber  and  a  Jones  matrix  $2^  for  propagation  in  the  op¬ 
posite  direction.  The  principle  of  reciprocity  relates  the  complex  scat¬ 
tering  coefficients.  The  representation  of  the  matrices  S12  and 
in  terms  of  the  complex  scattering  coefficients  are  presented  in  Eqs.  (la) 
and  (lb).  Equation  (Ic)  shows  the  relation  between  the  complex  scattering 
coefficients  as  a  result  of  reciprocity. 


5 


12 


t(xl ,x2) 

t(xl ,y2)\ 

t(yl ,x2) 

t(yl,y2)/ 

t(x2,xl ) 

t(x2,yl)\ 

t(y2,xl ) 

t(y2,yl)/ 

-  1  - 

(la) 


(lb) 


t(k,£) 


t  (Jl,k) 


k,£  =  xl,x2,yl,y2 


(Ic) 


Equation  (Ic)  means  that  the  direction  in  which  a  path  is  transversed  does 
not  matter.  For  a  lossless  fiber,  the  matrices  ^nd  S^-j  are  unitary. 

This  means  that  the  sum  of  the  absolute  squares  of  the  elements  in  any  row 
or  column  is  equal  to  1.  Loss  can  be  introduced  by  multiplying  the  matrices 

“CtL 

5^2  ^nd  S21  by  e  ,  where  L  is  the  length  of  the  fiber  and  a  is 
the  loss  per  unit  length  of  the  fiber.  In  the  analysis  of  the  polarization 
properties  of  various  fiber  optic  rotation  sensors  (FORS),  it  is  assumed 
that  the  fibers  are  lossless.  The  Sagnac  phase  shift  modifies  the  matrices 
S.|2  S21  as  shown  in  Eqs.  (2a)  and  (2b). 
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’21 


t(xl ,x2) 

t(xl ,y2)^ 

t(yl ,x2) 

t(yl,y2)i 

t(x2,xl ) 

t{x2,yl)\ 

t(y2,xl ) 

t(y2,yl)/ 

(2a) 


(2b) 


0^  is  the  Sagnac  phase  shift. 

Although  the  matrices  S.|2  and  $21  can  be  used  as  given  in  Eq.  (2a) 
and  (2b),  it  will  simplify  further  analysis  if  an  explicit  form  is  given 
for  the  t(k,Ji)'s.  The  t(k,£)'s  are  complex  numbers  which  represent  the 
effect  of  a  channel  on  a  light  wave,  propagating  from  one  port  to  another 
port  in  the  system.  They  will  be  represented  in  polar  forms.  The  magnitude 
of  t(k,J?,)  is  the  amplitude  of  the  electric  field  of  the  light  wave,  at  port 
k  for  unit  electric  field  applied  at  port  1.  The  phase  of  t(k,£)  is  the 
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phase  shift  undergone  by  a  light  wave  traveling  the  channel  from  port  I  to 
port  K.  The  2x2  matrix  of  the  t(k,£)'s  is  unitary.  Using  this  fact, 
the  t(k,Jl)'s  can  be  written  in  the  following  polar  forms. 


t(xl,x2)  =  cos  9.  exp[j0(xl ,x2)]  (3a) 
t(xl,y2)  =  sin  9,  exp[j0(xl ,y2) ]  (3b) 
t(yl,x2)  =  -sin  9.  exp[j0(yl ,x2}]  -  (3c) 
t(yl,y2)  =  cos  0.  exp[j0(yl ,y2)]  (3d) 


In  Eqs.  (3a),  (3b),  (3c),  (3d),  the  parameter  0  describes  the  rotation 
of  the  input  state  of  polarization  by  the  fiber.  The  phase  shifts  0(k,Jl) 
must  now  be  determined  by  considering  a  length  L  of  optical  fiber  with  a 
single  scattering  point  located  at  a  distance  Z  along  the  fiber.  This 
model  is  shown  in  Fig.  2.  The  phase  shifts  0(k,£)  calculated  from  this 
model  are  given  in  Eq.  (4). 


0(xl,x2)  =  3^  L 

(4a) 

HyUy2)  =  3y  L 

(4b) 

0(xl,y2)  =  3y  L  +  A6Z 

(4c) 

0(yl,x2)  =  3^  L  -  A3Z 

(4d) 

A3  -  6^-6^ 

(4e) 

These  equations  can  be  put  in  a  more  symmetric  form  by  factoring  out  the 
term  (3^  +  By)L/2  .  This  results  in  the  following  equations  for  the  phase 
shifts: 


0(xl  ,x2) 

=  3  L  +  A6L/2 

(5a) 

0(yi ,yz) 

=  3  L  -  ASL/2 

(5b) 
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Z=0  f  Z  Z=L 

I 

I 


Figure  2.  Schematic  representation  of  a  single  mode  fiber. 
The  paths  taken  by  the  two  polarizations  X^Y  are  shown 
explicitly  by  XX’  and  YY’  with  propagation  constants  0 
and  .  Z  can  represent  a  point  in  the  fiber  were  the 

power  in  the  two  modes  is  interchanged. 


0(xl,y2)  =  3  L  -  A6(L/2 

-  z) 

(5c) 

0(yl»x2)  =  6  L  +  A3(L/2 

-  z) 

(5d) 

3  =  (3^  +  3^)/2 

(5e) 

combining  Eqs.  (2),  (3),  and  (4),  an  explicit  representation  for  the  matrices 
$12  and  $21  is  found  for  the  case  of  a  fiber  with  a  single  scattering  cen¬ 
ter. 


’12 


’21 


/cosQ 

sin0  e" 

COS0  e” 

/CCS9 

-sinQ  e 

is1n9 

cos9  e"' 

jA3(L/2-Z)’ 


(6a} 


(6b) 


In  general  a  real  fiber  has  many  scattering  points.  To  find  the  matrices 
$12  and  $21  in  this  case,  divide  the  fiber  into  sections,  each  section 
containing  one  scattering  point,  then  multiply  the  matrices  for  each  section 
together  to  get  the  matrix  for  the  fiber.  The  matrices  $12  and  $21  have 
the  same  general  form  as  the  matrices  in  Eqs.  (6a)  and  (6b).  This  general 
form  is 
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21 


30^/2  j0(L) 
e  e 


-j03/2 


/  cos9 

\-sin9 

cos9 

/cos9 

j0(L)/2 

-sin9  e 

\sin0 

cos9 

(7a) 


(7b) 
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In  Eqs.  (7a)  and  (7b)  0(L)  is  an  average  phase  shift  which  depends  on 

the  length  of  the  fibers;  C  and  0  are  phase  shifts  which  depend  on 
the  birefringence  and  number  of  scattering  centers  in  the  fiber.  The 
physical  interpretations  of  C  and  0  are  as  follows:  0  is  the  nonre¬ 
ciprocal  phase  shift  introduced  by  the  fiber  and  its  value  is  the  phase 
difference  between  light  waves  which  travel  from  port  xl  to  y2  and  from 
port  x2  to  port  yl  (see  Fig.  1  of  Appendix  C);  ^  is  the  measure  of  the 
birefringence  of  the  fiber  and  its  value  is  the  phase  difference  between 
light  waves  which  travel  from  port  xl  to  port  x2  and  from  port  yl  to  y2. 

In  order  to  fully  represent  a  fiber  gyro,  it  is  necessary  to  develop 
Jones  matrices  for  polarizers  and  beam  splitters.  The  Jones  matrix  P  , 
for  a  general  polarizer  at  an  angle  a  to  the  x-axis  is  given  in  Eq.  (8). 

2  2 

p^cos  a  +  PySin  a  ^Px'Py^  (cosa) 

2  2 

(P^-Py)  (sin  a)  (cosa)  p^sin  a  +  p^cos  a 

In  this  expression,  p^  and  p^  are  the  fractions  of  the  x,y  electric 

fields  passed  by  the  polarizer.  For  an  ideal  polarizer  p  =  1  and  p  = 

A  y 

0  .  The  basis  polarizations  for  this  representation  of  a  polarizer  are 
the  X  and  y  linear  polarizations. 

A  beam  splitter  is  characterized  by  two  Jones  matrices,  one  for  trans¬ 
mission  (T)  and  one  for  reflection  (R)  .  The  form  of  these  matrices 

is  dependent  on  the  particular  basis  polarizations  chosen.  For  an  ideal, 
lossless,  50-50  beam  splitter  in  the  x  and  y  linear  polarization  basis, 
the  T  and  R  matrices  are  given  in  Eqs.  (9a)  and  (9b). 


T 


\  °  A 


(9a) 
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R 


'A  « 


(9b) 


The  T  and  R  matrices  which  represent  a  real  beamsplitter  in  the  x  and 
y  linear  polarization  basis  are  given  in  Eqs.  (10a)  and  (10b). 


(10b) 

In  Eqs.  (10a)  and  (10b),  t^  ,  t^  ,  r^  and  r^  are  complex  numbers  with 
1^2=1  and  =  l  . 

With  the  Jones  matrices  and  the  coherence  matrix  formalism  (Appendix  A) 
it  is  possible  to  analyze  various  types  of  fiber  gyros.  The  birefringence 
inherent  in  all  single  mode  optical  fibers  can  now  be  described  by  three 
measurable  quantities  0,0,  and  ^  .  This  can  result  in  better  charac¬ 
terization  of  fiber  birefringence  and  the  effects  of  fiber  birefringence  on 
fiber  gyroscopes. 
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There  appear  to  he  limitations  in  the  operation  t>f  optical-fiber  Sagnac  gyro  rotation  sensors  that  have  imposed  a 
minimum  measurable  rotation  rate  that  is  much  higher  than  that  caused  by  quantum  noise.  We  show  that  one 
source  of  limitation,  namely,  the  superposition  of  a  nonreciprocal  pair  of  waves  generated  by  backward  scattering 
from  the  incident  waves,  can  result  in  significant  error  hut  can  be  mitigated  by  appropriate  system  design  and  sig¬ 
nal  modulation. 


Fiber  gyro  progress  is  now  on  a  plateau  at  which  the 
observed  sensitivity  is  far  below  the  expected  range. 
Calculations  of  the  sensitivity  of  the  system,  assuming 
no  noise  sources  other  than  basic  quantum  noise,  indi¬ 
cate  a  capability  for  performance  extending  to  less  than 
10"'*  degree  of  rotation  per  hour.'  In  spite  of  the  at¬ 
tractive  promise,  no  experimental  fiber  gyro  has  come 
close  to  this  level.  Noise,  in  the  range  of  10^  to  10®  times 
higher  than  quantum  noise,  is  found  to  be  present.  The 
noise  is  sufficient  to  limit  the  sensitivity  to  the  range  of 
1.0  to  100  deg/h,  even  with  integration  times  extending 
to  many  seconds.-  ‘  Existing  analyses  of  the  system, 
covering  all  presumed  sources  of  noise,  indicate  no  un¬ 
avoidable  sources  that  have  a  magnitude  sufficient  to 
change  the  predicted  sensitivity.  In  this  Letter,  we 
identify  an  effect,  the  consequences  of  which  have  not 
been  treated,  give  estimates  of  its  magnitude,  and  pro¬ 
pose  new  practical  approaches  to  its  solution. 

Reciprocal  and  Nonreciprocal  Effects 

The  factors  that  previously  have  been  recognized  as 
limiting  the  accuracy  of  rotation  measurement  are, 
broadly,  residual  nonreciprocal  characteristics  in  the 
light  path,  nonhomogeneous  environmental  changes, 
quantum  noise  in  the  conversion  of  light  to  electrical 
signals,  and  electrical  noise  in  the  measuring  system. 
The  requirement  of  optical  reciprocity  in  the  two  waves 
has  dictated  the  use  of  single-mode  fiber,  in  order  to 
en.sure  that  there  are  identical  path.s  in  the  two  direc¬ 
tions,®'^  and  precise  environmental  control  of  inhomo¬ 
geneous  vibration  and  temperature  changes.*®  Quan¬ 
tum  noise  and  electrical  noise  have  been  considered,'* 
and  there  appear  to  be  adequate  sources,  fibers,  and 
detectors  for  the  precise  measurement  of  phase. 
Scattered  radiation  that  is  due  to  inhomogeneities  in  the 
fiber  produces  a  secondary  pair  of  waves  propagating 
in  the  fiber.  These  secondary  waves  and  the  mecha¬ 
nisms  producing  them  are  rather  well  known.  However, 
we  find  that  in  the  Sagnac  interferometer  these  same 
waves  can  produce  interfering  components  in  the  output 
that,  in  the  absence  of  conditions  that  randomize  their 
phase,  could  have  an  order  of  magnitude  consistent  with 
existing  sensitivity  limitations.  The  backscattered 
waves,  as  viewed  at  the  output  of  the  loop,  do  not  share 
the  property  of  reciprocity  possessed  by  the  original 


incident  waves  and  can  contain  a  component  that  can¬ 
not  be  treated  merely  as  white  noise. 

At  least  the  parts  of  the  scattered  waves  that  are  due 
to  physical  defects  in  the  fiber  and  to  Rayleigh  scat¬ 
tering  are,  to  some  degree,  temporally  and  spatially 
coherent  with  the  two  primary  waves  and  add  vecto- 
rially  at  unpredictable  phases  that  are  not  identical  for 
the  two  waves  (Fig.  1).  If  the  source  and  the  system 
were  absolutely  stable,  the  error  that  was  due  to  the 
secondary  waves  would  be  constant  and  could  be  cali¬ 
brated  without  limiting  the  sensitivity  or  accuracy  of  the 
system.  Conversely,  if  conditions  were  sufficiently 
variable,  coherence  could  be  completely  destroyed,  and 
the  waves  would  simply  contribute  shot  noise,  as  has 
previously  been  assumed.  In  fact,  however,  if  left  alone, 
the  secondary  waves  could  produce  an  unpredictable 
error  much  larger  than  shot-noise  contributions. 

The  forward  scattering  appears  to  have  little  effect 
on  the  uncertainty  in  phase  and  is  believed  to  be  re¬ 
ciprocal,  affecting  the  oppositely  directed  waves  alike. 
The  part  of  the  scattered  radiation  that  falls  within  the 
acceptance  angle  of  the  fiber  in  the  reverse  direction, 
however,  provides  the  secondary  waves  that  interfere 
with  the  desired  (primary)  counterpropagating 
waves. 

Quantitative  Evaluation 

An  estimate  of  the  magnitude  of  phase  uncertainty  in 
a  fiber  loop  can  be  made  from  simple  calculations.  As 
an  example,  consider  an  attenuation  exp(— a,L)  that  is 
due  to  Rayleigh  scattering  with  the  lost  energy  uni¬ 
formly  scattered  over  the  length  of  a  fiber  with  a  scat¬ 
tering  directivity  G  along  the  fiber.  The  fraction  of  the 
total  scattered  energy  that  is  captured  by  the  fiber  core 
is  where*^  /3  is  the  acceptance  angle  of  the  fiber 

core  and  I  <  G  <1.5,  The  detailed  values  of  G  and 
other  parameters  in  this  Letter  may  depart  from  the 
values  assumed,  depending  on  various  factors,  but  this 
departure  will  not  change  the  order  of  magnitude  of  our 
conclusions. 

For  a  loop  length  L,  it  can  easily  be  shown  that,  for 
small  attenuation,  the  ratio  of  the  power  in  the  back- 
scattered  secondary  wave  to  that  in  the  primary  wave 
after  one  circulation  is  given  by 
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SECONDARr  wave 

Fig.  1.  Fiber-loop  Sagnac  interferometer.  The  scatter  loss 
in  a  fiber  that  is  due  to  fixed  scattering  centers,  including 
Rayleigh  scattering  (the  principal  loss  mechanism  in  very- 
low-loss  fibers),  gives  rise  to  counterrotating  secondary  waves 
that  are  indistinguishable  from  the  signal  component  that  is 
caused  by  rotation.  Their  cumulative  effect  can  he  reduced 
by  suitably  modulating  the  input  signal.  <(>,  phase  change  that 
is  due  to  rotation  rate.  <t>',  phase  change  including  the  sec¬ 
ondary  wave  vectors. 


Table  t.  Parameters  for  Numerical  Estimate 


Symbols 

Parameters 

Values  Used 

X 

Wavelength 

1.0  nm 

c 

Free-space  velocity  of 
light 

3  X  10“^  msec 

n 

Refractive  index  of  fiber 

1.5 

N 

Nurnber  of  turns  in 
length  L 

318 

D 

Loop  diameter 

1-0  m 

Capture  angle  of  fiber 

0.1  rad 

G 

Directional  gain 

1.0 

Attenuation  constant 

1.2  X  lO-'m-i 
(~0.5  dB/km) 

L 

Fiber  length 

1000  m 

P,/Pi  (1) 

The  two  counterrotating  pairs  of  waves  can  combine 
to  give  a  maximum  phase  deviation  of 

=  2(P,/Pi)i/2  =  /3(Ga,L)i/2,  (2) 


which,  for  any  single  measurement,  is  indistinguishable 
from  a  phase  shift  induced  by  rotation.  The  rotation 
rate  that  corresponds  to  this  phase  deviation  is  obtained 
using  the  well-known  Sagnac  formula  and  is  given  by 


max 


Ac 

2(V(7rD)2 


(3) 


or,  substituting  N  =  L/(ttD), 


__  \c0  /G^,U/^ 

2TrD  L  I 


(4) 


Using  the  values  given  in  Table  1  for  a  single  transit 
around  a  multiturn  loop,  we  calculate  a  maximum  phase 
error  of  0.035  rad  that  linearly  translates  to  a  rotation 
error  of  up  to  0.095  deg/sec  or  341  deg/h. 


Minimization  of  Error  Due  to  Scattering 

There  are  a  number  of  ways  in  which  the  coherence 
between  the  primary  and  secondary  (scattered)  waves 
may  be  reduced.  Some  ways  may  reduce  the  reciprocity 


of  the  two  primary  waves,  but  some,  such  as  frequency 
modulation  of  the  primary  signal  or  a  physical  modu¬ 
lation  of  the  loop  length  in  a  controlled  fashion,  may  not. 
If  we  assume  such  a  modulation  at  a  level  high  enough 
so  that  the  number  of  wavelengths  in  the  loop  is  changed 
significantly  in  the  sampling  time,  the  secondary  waves 
add  to  the  emerging  primary  waves  with  random  phase. 
If  we  sample  at  a  frequency  of  M  per  unit  of  time,  and 
if  the  secondary  wave  phase  is  random  between  samples, 
the  uncertainty  in  the  angular  position  of  the  interfer¬ 
ometer  after  a  given  time  is  determined  by  a  random- 
walk  process  and  is  given  by 

e  =  ^  (5) 

In  the  foregoing  example,  taking  10  samples/sec  and 
integrating  for  a  full  hour  would  yield  a  random-walk 
rotation  error  of  1.27  deg/h^^^  This  error  is  of  the  same 
order  of  magnitude  as  that  of  the  best  operation  yet 
reported. 

The  error  that  is  due  to  scattering  is  probably  less¬ 
ened  somewhat  in  existing  fiber  gyros  because  of  certain 
factors  (some  inadvertent,  some  unavoidable)  that  act 
to  reduce  the  coherence  between  the  primary  and  sec¬ 
ondary  waves.  The  systems  referred  to  above  use 
modulation,  which  may  serve  to  randomize  the  phase 
of  the  scattered  waves  to  some  extent,  although  the 
modulation  is  employed  for  a  basically  different  pur¬ 
pose.  Some  degree  of  randomness  undoubtedly  occurs 
because  of  mechanical  and  thermal  changes,  hut  the 
changes  produce  other  errors  if  they  are  not  completely 
and  instantaneously  reciprocal.  Variation  in  laser 
frequency  is  another  source  of  randomization. 

These  factors  notwithstanding,  the  contribution  of 
backscattering  to  the  overall  error  may  still  be  signifi¬ 
cant  or  even  dominant.  Furthermore,  as  the  sensitivity 
of  these  systems  to  mechanical  and  thermal  perturba¬ 
tions  is  reduced,  an  increase  in  the  coherence  of  the 
scattered  waves  is  expected,  since  the  effects  of  the  re¬ 
maining  environmental  variations  are  in  general  too 
slow  to  permit  independent  random  sampling  at 
frequencies  high  enough  to  reduce  significantly  the  error 
that  is  due  to  scattering-  How  rapidly  can  we  sample 
and  expect  to  have  the  secondary  (scattered)  waves 
uncorrelated  in  phase?  It  will  take  additional  analysis 
to  determine  the  limits,  but  it  seems  probable  that  the 
effects  of  the  secondary  waves  can  be  made  very  small. 
For  instance,  if  we  use  a  pulsed  system,  with  pulses  in¬ 
troduced  at  a  rate  c/nL  (i.e.,  pulse  length  =  transit 
time),  and  shift  the  pulses  in  frequency  by  c/nL  per 
period  [which  is  (cfnL)-  per  second]  to  randomize  the 
phase  of  the  scattered  wave,  then  M  =  2  X  lO'"’  per  sec¬ 
ond.  The  error  that  is  due  to  Rayleigh  scattering  alone 
becomes  0.00015  deg/sec*''-  (1-sec  random  walk)  or  0.009 
deg/h^^“  (1-h  random  walk).  For  the  1-km  fiber,  this 
requires  a  200-kHz  per  pass  (40-GHz/sec  or  1.3-A/sec) 
change  in  source  frequency. 

The  ultimate  solution  for  the  minimization  of  this 
source  of  error  involves  reducing  the  magnitude  of  the 
backscattering  itself.  A  large  reduction  in  the  ampli¬ 
tudes  of  the  Rayleigh-scattered  waves  can  be  made  by 
using  short  pulses  (perhaps  1  nsec)  and  gating  the  out¬ 
put  so  that  only  the  scattered  energy  from  a  short 
(0.2-m)  length  of  fiber  contributes  to  the  part  of  the 
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secondary  wave  that  interferes.  The  short  length  that 
matters  could  be  specially  processed  or  constructed  to 
reduce  scattering.  This,  and  the  equivalent  use  of 
frequency  modulation  and  chirp  modulation,  are  the 
subjects  of  a  forthcoming  paper.  We  believe  that  im¬ 
provements  in  the  fiber  system  and  in  the  modulation 
proce.ss  can  reduce  this  source  of  error  to  a  level  less  than 
that  which  is  due  to  photon  noise  and  other  factors. 

Discussions  with  H.  Marcatili,  M.  Chodorow,  G. 
Pavlath,  and  H.  Lefevre  have  been  helpful  in  developing 
these  concepts. 

This  work  was  supported  by  the  U.S.  Air  Force  Office 
of  Scientific  Research  under  contract  F49620-80-C- 
0040. 
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